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Examples  Manual  for  Program  MAVAEIT 

E.  L.  Skiba,  G.  W.  McMahon  and  Z.  Wozniak 

ABSTEIACT 

This  document  is  the  third  in  a  set  of  documents  that  provides  information  on 
the  Finite  Element  Model  to  Analyze  the  Vibrations  and  Acoustic  Radiation  of 
Transducers  (MAVART).  The  other  two  documents  in  the  set  are: 

1 .  User's  Manual  for  Program  MAVART.  2.  Theoretical  Manual  for 
Program  MAVART. 

The  program  MAVART  has  been  developed  for  the  Defence  Research 
Establishment  Atlantic  (DREA)  tmder  research  contracts  to  Canadian  industry 
from  1976  to  the  present. 

This  manual  provides  a  number  of  sample  problems  to  aid  users  in  model 
development  and  output  interpretation.  It  demonstrates  virtually  all  elements 
and  features  of  MAVART.  It  also  provides  a  set  of  problems  for  validating 
modifications  and  new  installations. 

RESUME 

Ce  manuel  est  le  troisi^me  d'une  s6rie  de  documents  qui  contiennent  des 
renseignements  sur  le  logiciel  du  Module  d' analyse  des  rayonnements 
acoustiques  et  des  vibrations  des  transducteurs  par  616ments  finis  (MAVART). 
Les  de\ax  autres  documents  de  la  s6rie  sont:  1 .  le  Manuel  de  I'utilisateur  du 
logioiel  MAVART.  et  2.  le  Manuel  de  th4orie  du  logiciel  MAVART. 

MAVART  est  un  logiciel  qui  a  6t6  d6velopp6  pour  le  Centre  de  recherches  pour 
la  d^fense/Atlantique  dans  le  cadre  de  plusieurs  contrats  de  recherche  conclus 
avec  I'industrie  canadienne  de  1976  ^  aujourd'hui. 

Ce  manuel  donne  des  exemples  de  problfemes  congus  pour  aider  les 
utilisateurs  en  vue  du  developpement  de  modules,  et  de  I'interpr^tation  des 
r6sultats.  Ce  manuel  d4montre  pratisuement  tous  les  diverses  composantes  et 
fonctions  de  MAVART.  II  contient  aussi  une  s6rie  de  problfemes  pour  la 
validation  des  modifications  et  des  nouvelles  installations. 
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1  INTRODUCTION 

MAVART  is  a  two  dimensional,  axisymmetric,  dynamic  finite  element  code 
developed  for  Defence  Research  Establishment  Atlantic  (DREA)  tmder  research 
contracts  with  Canadian  industry.  Its  purpose  is  to  provide  a  computer  aided 
design  facility  for  axisymmetric  electroacoustic  transducers.  MAV ART's 
name  is  an  acronym  for  the  code's  function,  that  is,  a  Model  to  Analyze  the 
Vibrations  and  Acoustic  Radiation  of  Transducers. 

MAVART  was  first  implemented  in  1976,  under  a  contract  with  Acres 
International  Limited  of  Niagara  Falls,  Ontario.  It  has  been  revised  and 
improved  a  number  of  times  over  the  years,  both  at  Acres  and  elsewhere.  The 
current  version  (12.0)  resulted  from  a  contract  with  Acres  International. 

This  document  is  a  revised  edition  of  the  1987  Examples  Manual,  issued  as  part 
of  DREA  Contractor  Report  CR/87/442  [Ref.  1.1].  Companion  Theoretical  [1.1b] 
and  User's  Manuals  [1.1a]  complete  the  manual  set. 

The  purposes  of  this  revised  Examples  Manual  are: 

•  to  aid  MAVART  users  in  model  development,  and  data  input  and  output 
interpretation; 

•  to  provide  comprehensive  new  user  training; 

•  to  exercise  and  demonstrate  all  elements  and  features  of  MAVART; 

•  to  provide  a  set  of  problems  to  validate  proper  operation  of  the  program  after 
modifications  are  made; 

•  to  help  in  establishing  timing  benchmarks  for  reference  after  modifications 
are  made. 

In  addition  to  the  tests  reported  in  this  Examples  Manual,  some  special  testing 
of  MAVART's  capabilities  and  limitations  have  been  conducted.  These  tests 
are  reported  separately  and  are  referenced  in  this  document  when  they  are 
relevant  to  an  example  problem. 

All  of  the  programs  in  the  MAVART  package  are  written  in  Fortran  77,  and 
will  rim  on  DEC  VAX,  HP  9000,  Apple  Macintosh,  Silicon  Graphics,  and  IBM 
386  PC  compatible  platforms  without  further  modifications. 


1-1 


P153521.PDF  [Page:  16  of  124] 


P153521.PDF  [Page:  17  of  124] 


2  TEST  PROBLEM  DEFEMmON 

Most  test  problems  come  from  DREA  experience  with  simple,  typical 
transducer  designs  that  have  special  features  or  are  of  special  interest.  They 
have  been  selected  so  as  to  exercise  all  MAVART  features. 

2.1  LIST  OF  PROBLEMS 

The  problems  that  have  been  included  in  the  manual  are  grouped  below  by 
principal  element  type  and  are  individually  numbered. 

Series  1  Problems  PAR  Piezoelectric  Spherical  Shell 

1  STATC  pressure  load  (hydrophone  sensitivity,  stresses) 

2  EIGEN  resonance 

3  EIGEN  antiresonance 

4  CAP  AC  frequency  sweep 

5  CAP  AC  transient 

6  DRIVE  frequency  sweep  (stresses  at  resonance) 

Series  2  Problems  PT  Piezoelectric  Ring 

7  EIGEN  resonance  (  m  =  2,0  ) 

8  CAP  AC  frequency  sweep  (  m  =  0  ) 

9  DRIVE  resonance  search,  with  and  without  extra  DOF's 

carried 

Series  3  Problems  PAR,  SOLID  and  SLIDER  Trilaminar  Bender  Disk 

10  STATC  electric  stress 

11  STATC  thermal  stress 

12  CAP  AC  deformations,  with  and  without  slider  softening 

Series  4  Problems  SHELL  Thin  Disk/Tube 

13  EIGEN  flexural  and  extensional  resonances 

14  STATC  transverse  temperature  differential  and  nodal  load 

Series  5  Problem  Analytical  FLUID  and  FTOS 

15  DRIVE  piston  in  infinite  hard  baffle 

Series  6  Problem  SOLID  (CompKant  Fluid),  FLUID,  and  FTOS 

16  DRIVE  Compliant  fluid  modelling  of  a  resonating 

spherical  air  bubble 

Series  7  Problem  RING  and  SHELL 

17  EIGEN  RING-stiffened  SHELL  with  added  nodal  mass 
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Series  8  Problem  DAMPING  AppUcations 


18  DRIVE 
Series  9  Problems 

19  DRIVE 

20  DRIVE 

Series  10  Problem 


complex  modulus  and  viscous  damping  in  PT  ring 
Rigid  Node  Application 

Two  Free-flooding  PT  Rings  and  SOLID  Spacer 
RIGID  node  replacing  spacer 

Torsional  Applications 


21  STATC  SOLID  cylinder  with  torsional  loading  (m=0) 

22  EIGEN  SOLID  cylinder  in  torsional  resonance 


2.2  PROBLEM  /  FEATURE  MATRIX 

Figure  2.2.1  provides  a  cross  reference  chart  of  significant  program  features 
and  test  problems.  It  allows  the  MAVART  user  to  quickly  find  which  example 
problem  demonstrates  a  given  problem  t3q)e,  element  type,  node  type,  or  other 
program  feature.  The  listed  problems  attempt  to  exercise  all  of  the  features  to 
an  adequate  degree,  biased  towards  the  normal  use  of  the  program. 

Most  features  of  GRAFl  are  also  exercised  by  these  example  problems, 
although  no  cross  reference  chart  is  presented. 
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Figure  2.2.1:  Problem  /  Feature  Matrix 

FEATURE 

Problem  Niunber 

1  2  3  4  5  6  7  8  9  10  U12  13  14  15  16  17  18  19  2021222324 

p 

R  EIGE 
O  STAT 
B  CAPA 
I  DRIV 
D  FOUD 

Phase  TRAN 

E  SOLID(Q) 
L  PAR(Q) 

E  PT(Q) 

M  FLUID(Q) 
E  MEMB 
N  FTOS 
T  FTOF 
S  SILV 
SLIDER 
SHELL 
RING 

Vise 
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3  GENERAL  APPROACH  TO  PROBLEM  ANALYSIS 

The  problem  analysis  followed  the  sequence: 

•  preprocessing,  i.e.  finite  element  model  and  input  data  preparation; 

•  execution  in  MAVART;  and 

•  postprocessing,  to  interrogate  resxdts  of  the  analysis. 

For  this  edition  of  the  Examples  Manual,  the  problems  have  been  renm  on 
DREA's  |iVAX  3900.  Some  of  the  problems'  input  data  were  changed  slightly 
for  various  reasons,  and  these  changes  are  documented  where  they  occur. 

3.1  MODEL  DEVELOPMENT 

The  path  taken  to  prepare  the  finite  element  (F.E.)  model  for  a  problem  is: 

•  generation  of  the  finite  element  mesh  using  program  GRID; 

•  preparation  of  material  property  data  using  program  MATERL; 

•  addition  of  data  file  title,  control  and  other  data;  and 

•  minimization  of  the  model's  matrix  profile,  using  program  PROFIL, 
if  the  backing  store  solver  is  being  used  (not  required  for  SPARSPAKX 

In  addition,  the  mesh  may  be  modified  using  program  MOD,  for  subsequent 
parametric  studies. 


3.2  PROBLEM  SOLVING 

Each  MAVART  execution  produces  one  or  more  text  output  files  §RESxx.DAT, 
where  "§"  represents  a  user  identification  character  input  on  startup,  and  "xx" 
is  the  version  number.  A  binary  data  file  §PLTxx.DAT,  required  for  GRAFl 
input,  is  also  generated.  A  portion  of  the  data  in  §RESxx.DAT  is  output  to  the 
terminal,  or  to  a  Batch  LOG  file  so  that  the  operation  of  the  job  can  be 
monitored. 


3.3  DATA  PRESENTATION 

Each  MAVART  input  file  created  at  the  preprocessing  stage  bears  a  name 
compatible  with  the  problem  annotation,  e.g.  for  Problem  9  the  corresponding 
input  file  is  D09.DAT  .  Some  problems  have  multiple  versions  of  the  data  file, 
generally  with  letters  A,B,etc  appended  to  the  file  name. 

Each  section  of  the  following  chapters  contain  details  pertinent  to  a  particular 
problem  solved.  Significant  results  of  the  analysis  are  shown  in  the  form  of 
plots  generated  using  GRAFl,  such  as: 

•  model  grids; 

•  stress  and  strain  contours; 


1 


P153521.PDF  [Page:  22  of  124] 


•  pressure  and  pressure  gradient  contours; 

•  displacements  of  solid  boundaries; 

•  directional  response  of  acoustic  output;  and 

•  plots  of  various  parameters  versus  frequency. 

Most  of  the  graphic  output  presented  in  this  edition  was  generated  by  DREA's 
GRAFl  program,  but  some  minor  editing  was  done  in  the  MacDraw  II 
application  on  the  Macintosh  computer. 
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4  SERIES  1:  PAR  PIEZOELECTRIC  SPHERICAL  SHELL 

A  cross-sectional  diagram  of  the  radially  poled  piezoelectric  spherical  shell, 
made  of  Clevite  PZT-4  material  [4.1]  and  immersed  in  sea  water,  is  shown  in 
Figure  4.0.1.  The  inner  and  outer  radii  of  the  shell  are  0.911  m  and  1.0  m, 
respectively.  The  outer  radius  of  the  near-field  sea  water  sphere  used  in  the 
analysis  is  2.0  m. 

Due  to  symmetry  of  the  problem,  only  the  upper  hemisphere  has  been  analyzed. 
The  mesh  of  the  F.E.  model  is  shown  in  Figure  4.0.2  .  The  Z  and  R  axes  of  the 
model  system  of  cylindrical  coordinates  coincide  respectively  with  the  Z  and  X 
axes  of  the  coordinate  system  shown.  Element  numbers  are  printed  just  above 
the  centroid  of  the  element. 

Figures  4.0.3  (a,b)  and  4.0.4  (a,b)  give  details  of  node  and  element  numbering  of 
the  F.E.  model.  Note  that  the  fluid  node  numbers  are  printed  by  GRAFl  just 
above  the  node  position  and  solid  node  numbers  are  printed  just  below. 

The  elements  grouped  along  the  outer  boundary  of  the  fluid  are  FTOF  (Type  7) 
elements.  The  outer  surface  of  the  shell  is  "covered"  with  SILV  (T3rpe  8) 
elements.  These  are  necessary  for  antiresonance  analysis  of  a  PAR 
piezoelectric  structure,  but  can  be  useful  whenever  Type  A  (PAR)  nodes  must 
be  connected  together  electrically.  Also  along  the  outer  surface  of  the  shell  are 
FTOS  (T3T)e  6)  elements  that  connect  the  solid  structure  to  the  flxiid. 

The  horizontal  symmetry  plane  requires  proper  boundary  conditions  for  solid 
elements;  the  cross  plane  displacements  for  Nodes  1,3  and  4  are  restrained  in 
the  Z  direction,  i.e.,  JF]X(1)=1. 
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Figure  4.0.3  (a):  Series  1  Problems  -  Fluid  Node  Numbering 
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Figure  4.0.3  (b);  Series  1  Problems  -  Solid  and  Connecting  Fluid 
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Figure  4.0.4  (a):  Series  1  Problems  -  FLUID  and  FTOF  Element 
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Figure  4.0.4  (b):  Series  1  Problems  -  SOLID  and  SILV 
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4.1  PROBLEM  1 

Analysis  Type:  STATC  -  pressure  load  (hydrophone  sensitivity,  stresses) 

Input  File;  DOl  .DAT 

The  F.E.  model,  a  subset  of  the  full  model  from  Figure  4.0.1  ,  contains  PARQ, 
FTOS  and  SILV  elements  only.  The  nodes  used  are  shown  in  Figure  4.0.3(b). 

A  uniform  pressure  of  1  Pa,  applied  to  the  F  nodes  of  the  FTOS  elements,  is 
shown  in  input  data  as  a  pressure  DOF  fixity  and  shown  graphically  in  Figure 
4.1.1.  The  small  ''pins"  pointing  to  the  fixed  nodes  indicate  nonzero  fixities. 
Voltage  on  the  internal  surface  of  the  PARQ  elements  is  fixed  at  V  =  0.  GRAY 
and  SPIN  body  forces  (Card  13)  are  set  to  zero.  The  voltage  node,  to  which  the 
calculated  hydrophone  sensitivity  is  referred,  is  Node  5  and  input  pressure  is  1 
Pa  (Card  14).  It  is  important  to  notice  that  Node  5  is  located  on  the  silvered 
external  surface  of  the  piezoelectric  shell. 

Plots  of  selected  output  data  are  shown  in  Figures  4.1 .2  through  4.1 .4  .  Note 
that  the  radial  stress  depicted  in  Figure  4.1 .4  is  cylindrically  radial  and  not 
spherically  radial. 

The  hydrophone  receiving  sensitivity  printed  on  Figure  4.1.2  is  valid  at  low 
frequencies,  well  below  any  resonances.  Simple  theory  for  thin  spherical  shells 
leads  to  the  following  expression  for  receiving  sensitivity: 

M  =  rdzilE^  (4.1) 

where  r  =  0.9555  metre  is  the  mean  radius  of  the  spherical  shell, 

dzi  =  -123.E-12  mA^  is  the  transverse  piezoelectric  strain  constant, 

and  £33"^  =  1.15E— 8  F/m  is  the  dielectric  permittivity  at  constant  stress. 

With  these  "book"  values  for  PZT-4  [see  Ref.  4.1],  we  calculate  a  receiving 
sensitivity  of —159.98  dB  re  lV/|jPa,  in  good  agreement  with  the  value  obtained 
by  MAVART  of -160.00  dB. 
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Figure  4.1.3:  Problem  1  -  Axial  Stresses 
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Figure  4.1.4: 
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4.2  PROBLEM  2 

Analysis  Type;  EIGEN  -  resonance 

Input  FUe:  D02.DAT 

The  F.E.  model,  a  subset  of  the  full  model  from  Figure  4.0.1  ,  contains  PARQ, 
and  SILV  elements  only. 

As  the  resonance  condition  requires,  a  V  =  0  fixity  is  applied  to  all  of  the  nodes 
on  the  inner  and  outer  surfaces  of  the  shell  (i.e.,  a  short  circuit  condition).  Due 
to  the  presence  of  SILV  elements,  it  is  sufficient  to  apply  V  =  0  to  Node  4,  in 
order  to  assure  a  imiform  potential  of  V  =  0  on  the  entire  external  surface  of  the 
shell. 

Lower  and  upper  frequency  estimates  (Card  6)  are  set  to  300  Hz  and 
1000  Hz  respectively.  The  normalizing  DOF  number  was  chosen  as  83  after  a 
geometry  check  run  to  ensure  that  this  is  not  a  fixed  DOF.  The  value  of  the 
normalizing  DOF  is  set  to  l.OE— 6  .  Values  of  ERVAL  and  ERVEC  (Card  7)  were 
each  reduced  to  1  .OE— 6  to  ensure  convergence  to  an  accurate  solution  [See  Ref. 
4.2] 

Plots  of  mode  shape  for  the  first  four  modes  are  shown  in  Figures  4.2.1  through 
4.2.4.  Mode  4,  the  "breathing"  mode,  is  the  only  one  of  these  modes  that  can  be 
excited  piezoelectrically  with  the  (spherically)  radial  poling  condition  of 
Problem  2. 

Anal3dical  theory  [4.3]  for  a  thin  walled  spherical  shell  gives  the  following 
expression  for  the  resonance  frequency  (Mode  4): 

fr  =  r2  p  (sii  +  S12 

(4.2) 

where  r  is  the  mean  radius  of  the  sphere,  p  is  the  density,  and  sii  and  si2  are 
elastic  compliance  moduli  of  the  material.  The  short  circuit  elastic  moduli  are 
used  for  resonance  calculation,  and  open  circuit  moduli  for  antiresonance 
calculation,  giving  values  of  947  Hz  and  1162  Hz,  respectively,  using  Clevite 
"book  values"  for  the  properties  of  PZT-4  [4.1]. 
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Figure  4.2.1: 
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Figure  4.2.3:  Problein2  -  SrdMode  Shaiie 
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4.3  PROBLEMS 

Analysis  Type:  EIGEN  -  antiresonance 

Input  File:  DOS. DAT 

The  F.E.  model,  a  subset  of  the  full  model  from  Figure  4.0.1  ,  contains  PARQ, 
and  SILV  elements  only. 

For  antiresonance  analysis  an  open  circuit  condition  is  required;  hence,  a  V  =  0 
fixity  is  applied  to  one  siu-face,  in  this  case,  the  internal  surface  of  the  shell, 
while  the  outer  silvered  surface  is  free.  The  only  change  that  is  required  to 
Problem  2  input  data  is  to  remove  the  fixity  from  Node  4. 

The  modal  frequency  analysis  data  (Card  6)  are  the  same  as  for 
Problem  2.  Two  poling  conditions,  set  by  parameter  POLANG,  for  the 
piezoelectric  shell  have  been  used  in  both  resonance  and  antiresonance 
analyses:  (a)  Uniform  radial  (spherical)  poling  and  (b)  The  upper  five  of  the  ten 
PARQ  elements  have  their  poling  reversed.  The  first  four  mode  shapes  are 
essentially  identical  in  all  cases  to  those  shown  in  Figures  4.2.1  to  4.2.4  for  the 
resonance  (a)  condition.  The  frequencies  obtained  in  the  MAVART  analysis 
are  given  in  the  following  table  for  all  cases: 


Mode  # 

(a)  FULL  RADIAL  POLING 

Resonance  Antiresonance 

(b)  HALF  REVERSED  POLING 

Resonance  Antiresonance 

1 

404.37 

404.37 

404.43 

409.28 

2 

583.68 

583.68 

582.74 

583.02 

3 

838.12 

838.11 

844.18 

845.02 

4 

949.02 

1146.8 

948.54 

974.37 

Note  that  for  condition  (a),  the  antiresonance  frequency  differs  from  the 
resonance  only  for  Mode  4,  which  is  the  mode  of  most  interest  acoustically.  For 
condition  (b),  the  antiresonance  is  always  higher  in  frequency  than  the 
resonance,  indicating  that  the  other  modes  can  now  be  excited  piezoelectrically. 

The  resonance  and  antiresonance  frequencies  calculated  using  Eq.  4.2,  947  Hz 
and  1162  Hz,  respectively,  are  in  reasonable  agreement  with  the  MAVART 
predictions  of  949  Hz  and  1147  Hz,  above. 
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4.4  PROBLEM  4 

Analysis  Type:  CAP  AC  -  frequency  sweep 

Input  File;  D04.DAT 

The  F.E.  model,  a  subset  of  the  full  model  from  Figure  4.0.1,  contains  PARQ 
and  SILV  elements  only. 

A  V  =  0  fixity  is  applied  to  all  nodes  on  the  internal  surface  of  the  shell.  A 
sinusoidal  driving  voltage  of  1  V  amplitude  is  specified  as  a  fixity  at  Node  4  on 
the  silvered  surface. 

The  same  polarity  of  PARQ  elements  as  in  Problem  3,  condition  (b)  is  used,  so 
that  Mode  1  can  be  excited. 

A  frequency  scan  was  conducted  in  the  range  of  400  Hz  through  408  Hz,  with 
frequency  start  points  and  increments  specified  according  to  Card  8  format. 
Ten  element  groups  are  used,  one  element  per  group.  NSTAV,  the  'number  of 
staves'  in  the  transducer,  is  2,  as  a  hemisphere  is  being  analyzed;  for  a  full 
sphere,  NSTAV  would  be  1 . 

Plots  of  deformations,  obtained  at  404  Hz  excitation  frequency  for  two  different 
phase  angles,  are  shown  in  Figures  4.4.1  and  4.4.2.  Capacitance  is  plotted 
versus  frequency  in  Figure  4.4.3.  It  should  be  noted  that  the  capacitance  values 
passed  to  the  output  files  are  the  absolute  values.  It  is  necessary  to  use  absolute 
values  in  the  resonance  search  mode  (see  Problem  9)  since  the  logarithm  of  the 
value  is  used.  However,  any  negative  (i.  e.  inductance)  values  should  be 
restored  before  output.  This  will  be  corrected  in  a  future  revision  of  MAVART. 
Some  of  the  capacitance  values  just  above  the  resonance  in  Figure  4.4.3  may  be 
negative. 
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Figure  4.4.1: 
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4.5  PROBLEM  5 

Analysis  lype:  CAP  AC  -  transient 

Input  File:  D05.DAT 

Auxiliary 

Input  File:  D05FOUD.DAT 

Auxiliary 

Output  FHe:  FOUDC  .D  AT 

The  F.E.  model,  a  subset  of  the  full  model  from  Figure  4.0.1  ,  contains  PARQ 
and  SILV  elements  only. 

Voltage  fixities,  the  driving  point  and  the  polarity  of  PARQ  elements  are  the 
same  as  in  Problem  4  . 

In  order  to  obtain  a  time  dependent  solution  for  an  arbitrary  forcing  function,  it 
is  necessary  to  decompose  the  function  first  into  the  corresponding  Fourier 
series  using  the  Fourier  Decomposition  Preprocessor  invoked  by  setting 
variable  PROBED  to  "FOUD".  At  the  end  of  the  auxiliary  input  file 
D05FOUD.DAT,  the  Fourier  decomposition  data  are  given  in  a  format  specified 
by  Cards  16  (a,b)  .  In  this  case  there  are  sixteen  equally  spaced  values  in  one 
cycle  of  a  1-volt  amplitude  square  wave  (2  volts  peak-to-peak)  with  a 
fundamental  frequency  of  134.7  Hz.  Note  that  this  drive  is  not  a  true 
"transient",  but  a  repetitive,  nonsinusoidal  waveform. 

MAVART  execution  with  the  file  D05FOUD.DAT  as  input  yields  output  files 
§MAV.DAT  and  FOUDC.DAT.  The  latter,  containing  Fourier  coefficients,  is 
used  automatically  in  subsequent  MAVART  execution,  with  file  DOS  .DAT 
serving  as  the  main  source  of  input  data. 

In  file  DOS. DAT,  the  fundamental  frequency,  total  number  of  harmonic 
frequencies  and  number  of  times  at  which  displacements  and  stresses  are  to  be 
computed  are  given  according  to  Card  9  (a)  .  After  inspection  of  data  in  file 
FOUDC.DAT,  harmonic  numbers  1,  3  and  5  were  chosen  and  written  into 
D05.DAT,  as  required  by  Card  9  (b) .  Then,  six  equally  spaced  time  instances, 
ranging  from  7.5  ms  through  10  ms,  were  specified  as  per  Card  9  (c). 

Plots  representing  snapshots  of  the  deformation  at  the  selected  times  are  shown 
in  Figures  4.5.1  through  4.5.6  . 
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4.6  PROBLEM  6 

Analysis  Type:  DRIVE  -  frequency  sweep 

Input  File:  D06.DAT 

The  complete  F.E.  model  from  Figure  4.0.2  is  used. 

Voltage  fixities,  the  driving  point  and  the  polarity  of  PARQ  elements  are  the 
same  as  for  Problem  4  , 

The  frequency  sweep  is  conducted  in  only  one  range  of  325  Hz  through 

340  Hz,  with  a  rather  coarse  5  Hz  step.  Acoustic  radiation  loading  is  small  for 

this  mode  (Mode  1),  and  the  sweep  spans  the  narrow  resonance  region. 

The  output  data  indicate  that,  among  the  frequencies  surveyed,  330  Hz  is  closest 
to  resonance  in  sea  water.  Mass  loading  due  to  the  water  lowers  the  resonance 
frequency  from  the  in  air  value  of  404.4  Hz.  Figures  4.6.1  and  4.6.2  show 
frequency  dependence  of  admittance  and  transmitting  voltage  response, 
respectively.  A  cubic  spline  fit,  available  as  an  option  in  GRAFl,  has  been 
applied  to  these  data. 

Plots  of  selected  data,  obtained  for  330  Hz  excitation,  are  shown  in  Figures  4.6.3 
through  4.6.7  . 


4-19 


P153521.PDF  [Page:  43  of  124] 


Figure  4.6.3:  Problem  6  -  Displacement  of  Solid  Boundary 


PRC8IEM  #6,  PAR  Piezoelectric  Spherical  Shell 


FREQUENCY:  330.00  Hz 
•PROBLEM:  'DRIV 

•  RESISTANCE: 
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. REACTANCE : 
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Figure  4.6.4:  Problem  6  -  Directional  Response  at 330  Hz 


PROBLEM  #6, 
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PAR  Piezoelectric  Spherical  Shell 

with  reversed  polarity  on  half  of  shell  elements 
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Figure  4.6.5:  Problem  6  -  Near-Field  Pressure  Contours 


PROBLEM  #6,  PAR  Piezcelectric  Spherical  Shell;  reversed  poling  on  half 


ENTIRE  lEGION  RADIUS- 
ISOAMPLITDEE  CONTOURS 


2.000 


of  shell  elements 


FREQUENCY: 


330.00  Hz 
POSITIVE 
ZERO 
NEGATIVE 


ISOAMPLITUDE  CONTOURS 
.MAXIMUM  AMPLITUDE: 

77.66270  Pa 

.MINIMUM  AMPLITUDE: 

0.3  76  1  682  Pa 

.REFERENCE  AMPLITtDE: 


5.405019 

Pa 

AMPLITUEE  CONTOURS 
20.30000  dB 

17.40000 

dB 

14.50000 

dB 

11.60000 

dB 

8.700001 

dB 

5.800000 

dB 

2.900000 

dB 

O.OOOOOOOE+OO  dB 

-2.900000 

dB 

-5.800000 

dB 

-8.700001 

dB 

dB 

dB 

-17.40000 

dB 

-20.30000 

dB 

27 -OCT-8  9 

13:06 

Figure  4.6.6: 


Problem  6  -  Detail  of  Near-Field  Pressure 


PROBLEM  #6,  PAR  Piezcelectric  Spherical  Shell;  reversed  poling  cm  half 
Rl-  0.9000  R2-  1.5000  Zl-  0.7000  Z2=-  1.3000  of  shell  elements 


FREQUENCY:  330.00  Hz 
_  _  _  POgTIVE 

_ _ m  NEGATIVE 


isoampluude  ccntours 

.MAXIMUM  AMPLITUDE: 

11.83618  Pa 

.MINIMUM  AMPLITUDE: 

0.3  76  1  682  Pa 

.REFERENCE  AMPLITUDE: 

2.110070  Pa 

.AlyPLITUEE  CONTOURS: 
13.30000  dB 

11.40000  dB 

9.500000  dB 

7.600000  dB 

5.700000  dB 

3.8QOOOQ  dB 

1.900000  dB 

O.OOOOOOOE+OO  dB 
-1.900000  dB 

-3.800000  dB 

-5.700000  dB 

-7.600000  dB 

-9.500000  dB 

-11.40000  dB 

-13.30000  dB 


27-OCT-8  9  14:14 


ISOAMPLUUDE  CONTOURS 


4-22 


P153521.PDF  [Page:  45  of  124] 


Figure  4.6.7;  Problem  6  -  Near-Field  Pressure  Gradient  Contours 


PPOBIEM  #6,  PAR  Piezoelectric  Spherical  Shell;  reversed  poling  on  half 

of  shell  elements. 


STRESS  AND  STRAIN  PLOTTING 


PRESSURE  GRADIENT 

FREQUENCY:  340.00  Hz 

.DMENSICNS: 

(ENTIRE  REGION) 

RMCN-  0.0000000 
Rmx-  2.0000000 
ZMIN-  0.0000000 
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73.40000 
55.00000 
36.60000 
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5  SERIES  2:  PTPIEZOELECTmC  RING 

A  radial,  cross  sectional  diagram  of  the  piezoceramic  ring,  immersed  in  a 
sphere  of  sea  water,  is  given  in  Figure  5.0.1.  The  inner  radius  of  the  ring  is 
0.833  m,  the  outer  radius  is  1.0  m,  and  its  half  height  is  0.333  m.  It  is  made 
from  48  staves  of  tangentially  poled  Channel  5400.  The  outer  radius  of  the  sea 
water  sphere  is  2  m. 

The  mesh  of  the  F.E.  model  is  shown  in  Figure  5.0.2,  together  with  the  element 
numbers.  As  in  the  case  of  the  Series  1  model,  only  the  upper  hemisphere 
needs  to  be  modelled  because  of  S3mimetry.  Figures  5.0.3  (a,b)  give  details  of  the 
node  numbering  of  the  F.E.  model. 

The  ring  is  modelled  by  two  PTQ  elements.  The  fluid  model  contains  both 
quadrilateral  and  triangular  elements.  The  fluid  and  PT  solid  are  connected  by 
FTOS  elements.  FTOF  elements  are  grouped  along  the  outer  boundary  of  the 
fluid.  Node  303,  the  only  node  of  E  type,  is  positioned  at  the  origin  of  the  system 
of  coordinates  and  coincides  with  F  type  Node  236,  at  the  geometrical  centre  of 
the  model.  Coordinates  of  E  nodes  are  immaterial  in  the  analysis. 

The  horizontal  symmetry  plane  requires  proper  boundary  conditions  for  solid 
elements;  crossplane  (Z)  displacements  for  Nodes  98,  99  and  132  are  fixed  to 
zero. 
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Figure  6.0.3  (a):  Series  2  Problems  -  Outer  Node  Numbering 


Series  2  Problem^  PT  Piezoelectric  Ring 
GEOlffiTRY  PLOTTING 

NOCE  NtMBERS 

.DIMENSICKS: 

(ENTIRE  REGION) 
RMIN=  -0.1000000  m 
RMAX=  2.1000001  m 
ZMIW=  -0.1000000  m 
ZMAX=  2.1000001  m 

.MATERIAL  TYPES: 

1 

.ELEMENT  TYPES: 
4-6-7-14-15 
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Figure  5.0.3  (b):  Series  2  Problems  -  Inner  Node  Numbering  Detail 
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5.1  PROBLEM  7 

Analysis  Type;  EIGEN  -  resonance  (  m  =  2,0  ) 

Input  Files:  D07A.DAT  -  w  =  2 

D07B.DAT  -  m  =  0 

The  F.E,  model  is  a  subset  of  the  full  model  from  Figure  5.0.2  and  contains  the 
PTQ  elements  only. 

As  the  resonance  condition  requires,  a  V  =  0  fixity  is  applied  to  the  only 
electrical  DOF  at  Node  303  .  The  Fourier  mode  number  m  =  2  is  written  into  the 
problem  parameters  data  of  Card  2  (IMODE). 

Although  the  lower  and  upper  boxmds  on  the  frequency  estimate  were  set  to  100 
Hz  and  500  Hz  respectively,  the  predicted  resonance  frequency  was  found  at  72.9 
Hz,  outside  of  the  100  to  500  Hz  band.  Figure  5.1.1  shows  the  corresponding 
mode  shape  in  the  Z-R  plane  at  0  =  0.  Tangentially,  the  displacement 
dependence  is  proportional  to  cos  m0  (=  cos  20  in  this  case). 

The  same  input  data  file  was  run  with  Fourier  mode  number  m  set  to  zero, 
giving  a  resonance  frequency  of  516.06  Hz.  The  mode  shape  is  shown  in  Figure 
5.1.2. 

We  see  that  the  cross  sectional  area  of  the  ring  remains  essentially  constant 
under  vibration  when  m  =2,  indicating  the  pure  bending  nature  of  the  mode. 

On  the  other  hand,  when  m  =  0,  the  modes  are  extensional  and  the  area  varies 
in  vibration,  as  seen  in  Figure  5.1.2. 
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Figure  5.1.1:  Problem  7  -  1st  Mode  Shape  ( m  =  2 )  at  72.9  Hz 

ERCBLEM  #7,  PT  Piezoelectric  Ring  (m=2) 


DISPLACEMENT  CF  THE  SOLID  BOUNDARY 


Figure  5.1.2:  Problem  7  -  1st  Mode  Shape  ( m  =  0 )  at  516  Hz 
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5.2  PROBLEM  8 

Analysis  Typ©!  CAP  AC  -  frequency  sweep  (m  =  0) 

Input  File:  D08.DAT  (frequency  sweep) 

The  F.E.  model  is  a  subset  of  the  full  model  from  Figure  5.0.2  and  contains  the 
two  PTQ  elements  only. 

A  sinusoidal  driving  voltage  of  1  V  amplitude  is  specified  as  a  fixity  at  Node  303. 

Similar  to  Problem  4  ,  three  frequency  ranges  are  specified  for  scanning.  They 
cover  the  band  from  508  Hz  through  524  Hz.  One  element  group  composed  of 
both  PTQ  elements  is  used  (Card  11).  It  is  not  necessary  to  specify  the 
coordinates  of  an  integration  line  for  PT  elements;  i.  e.,  RRl,  ZZl,  RR2,  and  ZZ2 
are  not  used. 

In  the  output  data  the  maximum  computed  capacitance  was  found  at  a 
frequency  of  516  Hz.  As  this  is  very  close  to  the  actual  resonance  of  516.06  Hz 
predicted  in  Problem  7,  the  capacitance  is  very  high.  Being  an  undamped 
driven  problem,  the  capacitance  and  displacement  become  infinite  at  the 
resonance.  A  plot  of  the  deformation  at  516  Hz  is  shown  in  Figure  5.2.1  and  a 
plot  of  absolute  value  of  the  capacitance  versus  frequency  is  given  in  Figure 
5.2.2. 
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5.3  PROBLEM  9 

Analysis  Type:  DRIVE  -  resonance  search 

Input  Files:  D09A.DAT  -  extra  DOF  dropped 

D09B.DAT  -  extra  DOF  kept 

The  complete  F.E,  model  from  Figure  5.0.2  is  used. 

Two  versions  of  the  input  file  differ  in  the  extra  (w,  0  displacement)  degrees  of 
freedom  (DOF)  either  being  dropped  (IFOU=0),  or  retained  (IF0U=1),  in  the 
analysis.  The  extra  DOF's  are  carried  automatically  when  the  Fourier  mode 
number  m  >  0,  and  is  only  reqToired  for  m  =  0  if  data  from  different  Fourier 
modes  are  to  be  combined. 

The  voltage  fixity  is  the  same  as  in  Problem  8,  the  E  node  (#303)  is  fixed  at  1  volt. 

To  locate  the  resonance,  instead  of  frequency  sweep  as  in  Problem  2  ,  another 
strategy  is  applied  for  the  case  without  the  extra  DOF's.  Card  8  in  file 
D09A.DAT  specifies  only  one  frequency  range  within  which  the  search  is  to  be 
performed.  The  negative  value  of  parameter  F2  in  Card  8  activates  the 
resonance  search  mode,  which  continues  until  the  maximum  of  a  specified 
parameter  is  located  with  peak  tolerance  specified  by  F2I  in  decibels,  or  imtil  12 
iteration  cycles  are  reached.  If  the  maximum  is  at  one  or  the  other  end  of  the 
specified  search  range,  the  search  range  will  be  extended  until  a  lower  value  is 
seen  or  until  the  12  iteration  cycles  are  completed.  Various  parameters, 
including  a  DOF  at  a  specified  node,  may  be  chosen  as  the  search  parameter  by 
specif3dng  a  nonzero  value  for  NFUN  in  Card  12.  Allowable  values  for  NFUN  are 
discussed  in  Reference  [5.1],  In  a  DRIVE  analysis,  radiated  power  (i,  e., 
conductance,  NFUN  =  2)  is  the  default  parameter  on  which  the  resonance 
search  is  based;  in  CAPAC,  it  is  capacitance  (NFUN  =  3). 

The  parameter  NSEL  in  Card  12  specifies  a  node  whose  DOF  values  are  written 
to  the  plot  file  for  GRAFl  processing;  one  of  these  DOF's  may  be  used  in  the 
resonance  search  [5.1].  If  NSEL  is  zero,  the  default  node  used  is  the  first  in  the 
node  list. 

The  output  data  indicate  that  375  Hz  is  the  frequency  closest  to  resonance  in  sea 
water.  Plots  of  selected  output  data  are  shown  in  Figures  5.3.1  through  5,3.6. 
The  cubic  spline  fit  option  has  been  selected  in  GRAFl  to  plot  the  frequency 
response  data  in  Figures  5.3.3  and  5.3.4. 

The  case  with  the  extra  DOF's  carried  (D09B.DAT)  was  run  only  at  375  Hz 
excitation  frequency.  The  results  are  identical  to  the  previous  run  for  that 
frequency. 
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Figure  5.3.1:  Problem  9  •  Deformation  at 375  Hz 


PBOBLEM  #9,  PT  Piezoelectric  Ring  (m=0) 
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Figure  5.3.2:  Problem  9  -  Tangential  Stress  Contours 

PROBLEM  #9,  PT  Piezoelectric  Ring  (m=0) 


STRESS  AND  STRAIN  PLOTTING 
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Figure  S.3.5:  Pixjblein  9 -Directivity  at 375  Hz 
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Figure  5.3.6:  Problem  9  -  Near-Field  Pressure  Contours  at  375  Hz 
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6  SERIES  3:  TRILAMINAR  BENDER  DISK 

A  cross  sectional  diagram  of  the  trilaminar  disk  is  given  in  Figure  6.0.1.  It  is 
made  of  two  outer  layers  of  PZT-4  ceramic  1  cm  thick  with  a  2-cm  aluminiim 
layer  between.  The  piezoelectric  material  on  both  sides  of  the  disk  is  poled  in 
the  same  +Z  direction.  Although  the  structure  is  symmetrical  geometrically 
with  respect  to  the  R  axis,  it  cannot  be  considered  electrically  S3mimetrical  due 
to  the  chosen  piezoelectric  polarity,  and  as  a  result  the  complete  structure  has 
to  be  analyzed. 

SLIDER  elements  are  inserted  between  the  PZT-4  and  the  aluminum. 
Physically,  they  can  be  considered  as  representing  the  adhesive  joints  bonding 
the  materials  together,  and  are  only  necessary  if  the  adhesive  contributes  a 
significant  compliance  to  the  structure.  Otherwise,  the  model  can  be  simplified 
by  eliminating  the  SLIDERs  and  allowing  the  aluminum  to  share  the  adjoining 
A  nodes  of  the  PZT-4.  The  SLIDER  elements  can  serve  three  purposes: 

•  they  can  "soften"  the  junction  of  two  structural  components; 

•  they  provide  for  connection  of  two  dissimilar  nodes,  such  as  H  and  A  types, 
where  both  must  exist  because  of  the  elements  that  they  serve;  and 

•  they  can  be  used  to  apply  damping  in  complex  driven  analyses  (damping 
applications  are  covered  in  Problem  18) 

The  mesh  of  the  F.E.  model  is  shown  in  Figure  6.0.2.  There  was  no  necessity  of 
using  triangular  elements  to  model  the  structure;  they  were  chosen  in  order  to 
demonstrate  their  use. 

Figures  6.0.3  and  6.0.4  (a,b)  give  details  of  element  and  node  numbering, 
respectively.  The  axial  scale  has  been  expanded  to  allow  the  numbers  to  be  read 
more  easily. 


Figure  6.0.1:  Series  3  Problems  -  Cross-Sectional  Diagram 
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Figure  6.0.2:  Series  3  Problems  -  F.E.  Model 


Series  3  Problems,  PfiR  Trilaminar  Bender  Disk 


GECMETRY  PLOTTIN3 


ELEMENT  PLOTTING 

■  DIMENSICNS: 
(ENTIRE  REGION) 
RMIN=  -0.0034000 
RMRX=  0.2034000 

ZMIN - 0.1034000 

ZMAX”  0.1034000 

.MATERIAL  TYPES: 
2-3-9 

.ELEMENT  TYPES: 
1-2-9 


1  2  3  4  5  6  7  8  9  10  28-NOV-89 


Figure  6.0.3:  Series  3  Problems  -  Element  Numbering 


Series  3  Problems,  PAR  Trilaminar  Bender  Disk 


GEOMETRY  PLOTTING 


ELEMENT  NUMBERS 

.DIMENSICNS: 
(ENTIRE  REGION) 
RMIN=  -0.0034000 
RMAX=  0.2034000 

ZMIN - 0.0413600 

ZMAX=  0.0413600 

•MATERIAL  TYPES: 
2-3-9 

.ELEMENT  TYPES: 
1-2-9 


28-NDV-89  09:02 
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Figure  6.0.4(a):  Series  3  Problems  -  Node  Numbering,  Aluminum  Layer 

Series  3  Problems,  P®  Trilaminar  Bender  Disk 


GEOMETRir  PLOTTING 


Figure  6.0.4  (b):  Series  3  Problems  -  Node  Numbering,  Piezoelectric  Layers 

Series  3  Problems,  PSR  Trilaminar  Bender  Disk 


GEOMETRY  PLOTTING 
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6.1  PROBLEM  10 

Analysis  Type^  STATC  -  electric  stresses 

Input  File:  DIO.DAT 

The  complete  F.E.  model  from  Figure  6.0.2  is  used. 

A  1-volt  fixity  is  applied  to  all  A  nodes  on  the  top  and  the  bottom  surfaces  of  the 
disk.  A  zero  volt  fixity  is  applied  to  all  A  nodes  attached  to  the  SLIDER 
elements.  The  SLIDER  stiffness  is  set  at  l.OE+15  Pa  in  all  three  moduli.  Node 
77  is  restrained  in  the  Z  direction. 

A  plot  of  the  deformation  of  the  disk  is  shown  in  Figure  6.1.1.  In  Figure  6.1.2, 
an  expanded  axial  scale  is  used  to  show  radial  strain  contours  in  the  materials; 
the  solid  contour  line  is  the  line  of  zero  strain,  long  dashed  lines  are  positive 
strain  contours,  and  the  short  dashed  lines  are  negative  contours.  The 
waviness  in  the  contours  is  an  artifice  caused  by  the  use  of  the  triangular 
elements. 
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Figure  6.1.1:  Problem  10  -  Electrical  Deformation  with  a  Hard 

SIJDEtt 


PSR  Trilaminar  Bender  Disk,  Electrically  stressed, SLIDER  1.0E15  Pa 


DISPLACEMENT  CF  THE  SOLID  BOUNDARY 


Figure  6.1.2: 


Problem  10  -  Radial  Strains 


PEDBIEM  10  -  PAR  Trilaminar  Bender  Disk,  Electrically  stressed 

.STRESS  AM)  SIRAIN  PLOTTING  'STATC 

R  STRAIN 


V 


.DIMENSIONS: 

(ENTIRE  REGION) 

RMIN - 0.0034000  m 

RMAX=  0.2034000  m 
ZMIN=  -0.0440000  m 
ZMAX=  0.0440000  m 

.MATERIAL  TYPES: 
2-3-9 


.MAX  VALOE:0. 1340E-07 
AT  NODE:  171 
.VALUE  CONTOURS  AT: 
0.1200000E-07 
0. 80000QOE-08 
0.4000000E-08 
Q.OOOOOOOE+00 
-0.4000000E-08 
-0.8000000E-08 
-0.1200000E-07 


10  27-NCV-89 


13:  07 
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6.2  PROBLEM  11 

Analysis  Type:  STATC  -  thermal  stress 

Input  File:  Dll  .DAT 

The  complete  model  from  Figure  6.0.2  is  used. 

A  zero  volt  fixity  is  applied  to  all  A  nodes  on  the  top  and  bottom  surfaces  of  the 
piezoelectric  layers.  The  SLIDER  stiffness  is  set  at  l.OE+15  and  Node  77  is 
restrained  in  the  Z  direction,  as  in  Problem  10. 

The  element  data  contains  extension  Cards  5  (b)  specifying  element 
temperature  differences  from  zero  strain:  +20°  C  for  elements  of  the  top 
piezoelectric  layer,  —20°  C  for  the  bottom  piezoelectric  layer.  The  temperature  of 
the  aluminum  disk  is  left  unchanged  from  its  zero  strain  state. 

Figures  6.2.1  and  6.2.2  show  output  data  for  the  thermally  stressed  disk  that 
corresponds  to  Figures  6.1.1  and  6.1.2  for  the  electrically  stressed  disk.  We  note 
that  the  form  of  the  deformation  is  essentially  the  same  in  both  cases  but  the 
strain  is  very  different:  there  are  three  lines  of  zero  strain  (solid  lines)  with 
thermal  stresses  but  only  one  with  electrical.  Adjacent  parts  of  the  materials 
undergo  opposite  strains.  The  hydrophone  sensitivity  printed  on  Figure  6.2.1  is 
meaningless  in  this  situation. 
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Figure  6.2.1: 


Problem  11  -  Tbermal  Deformation  wiHi  a  Hard  SLIDER 


PBDBIEM  11  -  PflR  Trilaminar  Bmder  Disk,  thermal  stresses 
DISPLACEMENT  CF  THE  SOLID  BOUNDARY 


10- 


Figure  6.2.2: 


■g - T 


Problem  11  -  Radial  Strains 


.  PBDBIEM :  ’  STAT' 

H/P  SENSITIVITY: 
-542.02ciB  (volts/uPa) 

•  DISPIACEMENT  SCALE  R 

I  u:’JJB"7'y92E-U3  m — I 

.DISPLACEMENT  SCALE  Z 

l“D';'33r7-g9'2E-05' m  '  i 

RMIN=  -0.0019200  m 
RMAX=  0.206  0  800  m 
ZMIN=  -0.1Q19200  m 
ZMAX=  0.1060800  m 


10  27-Ni:V-89  11:46 


PROBLEM  11  -  PAR  TrJ.jaminar  Bender  Disk,  thermal  stresses 


R  STRAIN 

.DIMENSICNS: 

(ENTIRE  REGION) 
RMIN“  -0.0034000  m 
RMAX=  0.2034000  m 
ZMIN-  -0.04  40000  m 
ZMAX=  0.0440000  m 


.MATERIAL  TYPES: 
2-3-9 


•  MAX  VALDE:0.3957E-04 
AT  NODE:  18  9 
.VALUE  CCNTOURS  AT: 

0. 300000 QE -04 
0.2000QOOE-04 
0. 100000  OE -04 
0. 000000 OE +00 
-0.1000000E-04 
-0. 200000  OE-04 
-0.3000000E-04 


27-NOV-89  13:00 


STRESS  AND  STRAIN  PLOTTING 
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6.3  PROBLEM  12 
Analysis  Typ©!  C  APAC 


Input  Files: 


D12A.DAT 

D12B.DAT 

D12C.DAT 

D12D.DAT 


-  without  slider  softening  (a) 

-  with  full  slider  softening  (b) 

-  with  partial  slider  softening  (c) 

-  hard  slider 

-  soft  slider 

-  softened  slider 

-  hard  slider  resonance 


The  complete  F.E.  model  from  Figure  6.0.2  is  used. 


A  sinusoidal  100-Hz  drive  with  an  amplitude  of  1  volt  is  applied  to  all  A  nodes 
on  the  top  and  bottom  surfaces  of  the  disk.  All  A  nodes  attached  to  the  SLIDER 
elements  share  a  0  volt  fixity.  No  structural  nodes  are  restrained  in  the 
Z  direction. 


Two  element  groups,  one  for  each  piezoelectric  layer,  are  used.  Three  SLIDER 
conditions  are  analysed.  The  material  properties  of  the  SLIDER  elements  in 
Case  (a)  make  the  aluminium-piezoelectric  connection  stiff  (cn  =  C22  =  C33  = 
l.OE-i-15),  as  in  Problem  10.  In  Case  (b)  the  connection  remains  stiff  in  the  . 

Z  direction  while  allowing  free  sliding  in  the  R  direction  (022  =  C33  =  0.).  Case  (c) 
is  an  intermediate  condition  (c22  =  C33  =  l.OE+11),  providing  some  shear 
compliance  in  the  R-0  plane. 

Plots  of  deformations  at  100  Hz  are  shown  in  Figirres  6.3.1  to  6.3.3  for  the  three 
cases  (a),  (b),  and  (c),  respectively.  The  clamping  effect  of  the  aluminum  plate 
on  the  piezoelectric  layers  is  reflected  in  the  reduced  capacitance  as  the  sliders 
become  stiffer. 


An  'EIGE'  run  was  carried  out  for  the  hard  slider.  Case  (a),  predicting  a 
resonance  of  1647.8  Hz.  Deformations  are  shown  in  Figure  5.3.4.  Since  no  node 
is  fixed,  the  first  mode  will  always  be  at  zero  frequency,  hence  the  data  indicate 
a  NATURAL  MODE  NUMBER  of  '2'. 
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Figure  6.3.1:  Problem  12  -  Deformations  without  Slider  Softening 


PROBLEM  12  -  PfiR  Trilaminar  Bender  Disk,  hard  SLIDER,  l.OE+15  Pa 


DISPLACEMENT  CF  THE  SOLID  BOUNDARY 


■  DRIVEN  FREQOEICY: 

100.00  Hz 

■  PRDBIEM:  'CAPA' 

.  CAPACITANCE : 

0.274  9  600E-06  F 

■  DISPLACEMENT  SCALE  R 

I  0.2S7122QE-07  m  I 

■  DISPLACEMENT  SCAIE  Z 

I  U.2b/i22UE-U  I  m  1 

RMIN - 0.0040000  m 

RMAX=  0.22  7  4  00  0  m 
ZMIN“  -0.1040000  m 
ZMAX=  0.1274000  m 


l-DEX:-89  13:01 


Figure  6.3.2:  Problem  12  -  Deformations  with  Slider  Softening 


PROBLEM  12  —  PAR  Trilaminar  Bender  Disk,  soft  SLIDER 
DISPLACEMENT  CF  THE  SOLID  BOUNDARY 
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Figure  6.3.3:  Problem  12  >  Deformations  with  Partial  Slider 

Softening 

PK5BLEM  12  -  PAR  Trilaminar  Bender  Disk,  softened  SLIDER,  l.OE+11  Pa 

•DRIVEN  FREQUENCY: 
100.00  Hz 

•  PRDBIEM:  'CAPA' 

.  CAPACITANCE : 
0.2786057E-06  F 

•  DISPLACEMENT  SCALE  R 
I  0.2492C)6SS-(3'7  m  I 

•  DISPIACEMENT  SCALE  Z 

l"U.2  4y2UbJE-U  /  m  I 

RMIN-  -0.0040000  m 
RMAX-  0.2274000  m 
ZMIN=  -0. 1040000  m 
ZMAX-  0.1274000  m 


l-DBC-89  13:03 


Figure  6.3.4:  Problem  12  -  First  Mode  Shape 

PROBLEM  12  -  PAR  Trilaminar  Bender  Disk,  hard  SLIDER,  Resonance 


^DISPLACEMENT  (F  THE  SOLID  BOUNDARY 

5 

S 


3 

% - 1 - 5 - 5 - - S - E - 7 - § - § - r 
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7  SERIES  4:  SHELL  THIN  DISK)TUBE 

This  series  of  problems  demonstrates  the  application  of  SHELL  (Type  10) 
elements  in  modelling  a  thin  aluminum  disk  and  a  thin  walled  aluminum 
tube.  A  cross  sectional  diagram  of  the  disk  is  given  in  Figure  7.0.1  and  of  the 
tube  in  Figure  7.0.2.  The  radius  of  the  disk  is  10  cm,  and  the  tube  is  10  cm  long 
and  5  cm  in  radius.  In  both  cases  the  shell  thickness  is  2  mm,  and  is  manifest 
in  the  material  properties  as  discussed  in  the  Theoretical  Manual  [1.2].  The 
properties  specified  for  aluminum  in  program  MATER  (See  [1.6],  [1.7])  are 
Yormg’s  modulus  =  7.16E10  Pa,  density  =  2700  kg/m^,  Poisson's  ratio  =  0.344 
and  coefficient  of  thermal  expansion  =  23.4E— 6  /°C. 

The  F.E.  model  of  the  disk  with  node  and  element  numbering,  is  shown  in 
Figures  7.0.3  and  7.0.4,  respectively.  The  corresponding  F.E.  model  for  the  tube 
is  given  in  Figures  7.0.5  and  7.0.6.  GRAFl  has  no  provision  for  indicating  the 
thickness  of  a  shell. 
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Figure  7.0.3:  Series  4  Problems  -  Disk  Node  Numbering 


Series  4,  SHELL,  thin  disk 


GEOMETBY  PLOTTINS 


NOEE  NUMBERS 

•DIMENSICNS: 

(ENURE  REGION) 
RMIN=  -0.0020000  m 
RMAX=  0.1020000  m 
ZMIN=  -0.0520000  m 
ZMAX=  0.0520000  m 

•  MATERIAL  TYPES; 

10 

.ELEMENT  TYPES: 

10 


6-EEC-89  15:4  7 


Figure  7.0.4:  Series  4  Problems  -  Disk  Element  Numbering 


Series  4,  SHELL,  thin  disk 


GEOMETRY  PLOTTIIE 


ELEMENT  NUMBERS 

•  DIMENSICKS  : 

(ENTIRE  REGION) 
RMIN=  -0.0020000  m 
RNAX=  0.1020000  m 
ZMIN=  -0.0520000  m 
ZMAX=  0.0520000  m 

.MATERIAL  TYPES: 

10 

.EIEMENT  TYPES: 

10 
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Figure  7.0.5:  Series  4  Problems  -  Tube  Node  Numbering 


NDDE  NOMBERS 

DIMENSICNS: 

(ENTIRE  RESIGN) 
RMIN=  -0.0050000  m 
RMAX-  0.0830000  m 

ZMIN - 0.0550000  m 

0.0550000  m 

.MATERIAL  TYPES: 

10 

.EIEMENT  TYPES: 

10 


7-DB:-89  15:27 


Figure  7.0.6:  Series  4  Problems  -  Tube  Element  Numbering 


ELEMENT  NUtBERS 

.DIMENSICNS: 

(ENTIRE  REGION) 
RMIN=  -0.00500  00  m 
RMAX=  0.0830000  m 
ZMIN=  -0.05500  00  m 
ZMAX=  0.0550000  m 

.MATERIAL  TYPES: 

10 

.EIEMENT  TYPES: 

10 


7-DBC-89  15:26 


Series  4,  SHELL,  thin-walled  tube 
GEOMETRY  PLOTTING 


Series  4,  SHELL,  thin-walled  tube 
GEOMETRY  PLOTTING 
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7.1  PROBLEM  13 

Analysis  Type:  EIGEN  -  flexural  (disk)  and  hoop  (tube)  modes 

Input  FHes:  D13A.DAT  -  disk 

D13B.DAT  -  tube,  ends  restrained  (infinite  tube) 

D13C.DAT  -  tube,  centre  restrained  (finite  tube) 

D13D.DAT  -  tube,  hoop  mode 

All  entries  into  the  material  property  table  of  the  aluminum  SHELL  elements, 
according  to  the  data  requirements  [1.1a],  [1.7],  depend  on  the  shell  thickness. 

For  the  disk,  the  search  for  resonances  is  confined  to  the  frequency  band 
between  350  Hz  and  450  Hz.  As  in  Problem  2,  ERVAL  and  ERVEC  (Card  7) 
were  each  reduced  to  l.OE-6  to  ensure  accurate  convergence  [4.2].  The 
resonance  frequency  was  located  at  458.35  Hz.  A  plot  of  the  corresponding 
mode  shape  is  shown  in  Figure  7.1.1,  where  the  vibration  is  displayed  at  five 
phase  angles,  71/4  radians  apart.  (The  27c/4  vibration  curve  lies  along  the  disk,  so 
is  not  seen.)  Further  testing  of  flexural  and  extensional  modes  in  the  thin  disk 
are  reported  in  [4.2] 

For  the  tube,  since  hoop  type  modes  are  expected  to  fall  within  a  much  higher 
frequency  range,  a  15,000  Hz  through  16,500  Hz  band  is  used.  Also,  a  nrunber 
of  modes  are  expected  in  close  proximity,  so  the  tolerance  for  convergence  of 
eigenvalues  (EKVAL  on  Card  7)  is  lowered  to  1  .OE-8  and  the  tolerance  for 
convergence  of  eigenvectors,  ERVEC  is  l.OE— 6.  Similarly,  in  anticipation  of  an 
increase  in  the  number  of  necessary  iteration  cycles,  the  "stop  iteration 
number",  ITSP  is  increased  to  50. 

When  the  ends  of  the  tube  are  fixed  in  Z  and  p  (rotations),  two  modes  were 
found  in  the  search  range.  The  first  mode  is  antisymmetrical  at  15,839  Hz  and 
the  second  is  symmetrical  at  16,403  Hz.  A  plot  of  the  first  mode  shape  is  shown 
in  Figure  7.1.2.  Extending  the  search  range,  we  find  a  pure  hoop  mode  at 
17,456  Hz  as  Mode  4.  The  end  fixities  simulate  the  constraints  in  an  infinite 
tube,  and  calculation  of  this  resonance  from  simple  elastic  theory  gives  17,457 
Hz. 

When  there  are  no  fixities,  three  modes  occur  in  the  search  range,  two 
symmetrical  modes  at  15,838  and  16,396  Hz,  and  one  antisymmetrical  mode  at 
16,314  Hz.  The  first  mode  shape  is  shown  in  Figure  7.1.3. 

In  the  free  tube,  no  real  mode  displays  pure  hoop  vibration,  but  when  rotations 
are  suppressed  by  artificially  applying  |3  fixities  of  zero  to  all  nodes,  a  hoop 
mode  is  found  at  15,924  Hz.  The  mode  shape,  shown  in  Figure  7.1.4,  displays 
both  radial  and  axial  vibration.  If,  in  addition  to  the  (3  fixities,  the  axial  motion 
is  suppressed  by  Z  fixities  applied  to  the  ends  of  the  tube,  the  pure  hoop  mode  at 
17,456  is  again  observed,  now  as  Mode  1. 
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Figure  7.1.1;  Problem  13  -  Disk,  First  Flexural  Mode  Shape 

PROBLEM  #13A,  SHELL,  thin  disk 


DISPLACEMENT  CF  THE  SOLID  BOUNDARY 


Figure  7.1.2:  Problem  13  -  Tube,  First  Mode  Shape,  Fixed  Ends 


PROBLEM  13B,  SHELL,  thin-walled  tube,  ends  fixed  in  Z  and  p 


.NATURAL  FREQUENCY: 
15838.53  Hz 

.PRDBIEM:  'EIGE' 

.FOURIER  MCOE  NUMBER: 
0 

.NATURAL  MCOE  NUMBER: 
1 

.EIGEN  VAL  DEVIATION: 
0.  7501590E-09 

.NUMBER  OF  ITERATIONS 
33 

.DISPIACEMENT  SCALE  R 

I - 3.7S119-} - nH 

.DISPIACEMENT  SCALE  Z 

I 3.  /bll94 - m— I 

RMIN-  -0.0150000  m 
RMAX=  0.0890000  m 

ZMIN - 0.0520000  m 

ZMAX-  0.0520000  m 


6-D0C-89  16:10 


DISPIACEMENT  CF  THE  SOLID  BOUNDARY 


- 9 - Vi 
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Figure  7.1.3:  Problem  13  -  Tube,  First  Mode  Shape,  Free  Ends 


PROBLEM  13C,  SHELL,  thin-walled  tube,  centre  node  fixed  in  Z 


.NATHPAL  FREQUENCY: 
15838.00  Hz 

.PRDBIEM:  'EIGE' 

.FOURIER  MODE  NUMBER: 
0 

•NATURAL  MODE  NUMBER: 
1 

.EIGEN  VAL  DEVIATION: 
0.120  8  875E-09 

.NUMBER  OF  ITERATIONS 
3 


.DISPIACEMENT  SCAIE  R 


RMIN-  -0.0150000  m 
RMAX=  0.0890000  m 
ZMIN=  -0.0533000  m 
ZMAX-  0.0533000  m 


6-033-89  16:07 


Figure  7.1.4;  Problem  13  -  Tube,  Hoop  Mode  Shape,  Fixed  Rotations 


PRDBIEM  13D,  SHELL,  thin-walled  tube,  all  rotations  fixed 


NATURAL  FREQUENCY : 
15924.41  Hz 

.PROBLEM:  'EIGE' 

.FOURIER  MODE  NUMBER: 
0 

.NATURAL  MODE  NUMBER: 
2 

.EIGEN  VAL  DEVIATION: 
0. 583194  6E-07 

.NUMBER  OF  ITERATIONS 
8 


.DISPIACEMENT  SCALE  R 


.DISPIACEMENT  SCALE  Z 


RMIN^ - 0.0150000  m 

RMAX=  0.0890000  m 
ZMIN=  -0.0520000  m 
ZMAX-  0.0520000  m 


7-DEC-89  15:15 
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7.2  PROBLEM  14 

Analysis  Type:  STATC  -  transverse  temperature  differential  (a), 

-  and  nodal  load  (b) 

Input  Files;  P14A.DAT 

P14B.DAT 

Only  the  disk  structure  is  used  in  these  problems. 

In  Case  (a)  a  20°  C  transverse  temperature  difference  is  assumed  across  the 
shell  thickness,  with  the  higher  temperature  on  the  top  surface  of  the  disk.  The 
data  format  requires  that,  for  the  right  side  of  the  element  (as  seen  from  the 
first  node  on  the  element)  being  hot,  the  temperature  difference  is  positive.  In 
this  case  the  left  side  is  hotter,  so  a  negative  value  of  DELDT  is  used  on  the 
element  extension  data.  Card  5b. 

A  plot  of  the  thermally  deformed  disk  is  shown  in  Figure  7.2.1. 

In  Case  (b)  the  disk  is  clamped  around  its  periphery  by  specifying  zero  u,v,w, 
and  P  fixities  at  Node  21  .  A 10-N  force  is  applied  as  a  load  fixity  at  Node  1  in  the 
+Z  direction. 

A  plot  of  the  deformed  disk  is  shown  in  Figure  7.2.2.  In  both  Figures  7.2.1  and 
7.2.2,  the  printed  hydrophone  sensitivity  is  meaningless. 
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ttions 


.PFDBIEM:  'SIAT' 

H/P  SENSITIVITY: 
-194.27dE  (vclts/uPa) 

■  DISPIACEMENT  SCALE  R 

I  U.  ybBib-lUE-UJ  m'  I 


•  DISPLACEMENT  ^AIE  Z 

I  0.9581S40E-03  m  1 

RMIN=  -0.0020000  m 
RMAX=  0.1020  000  m 

ZMIN - 0.0520000  m 

ZMAX=  0.0520000  m 


b  8-DEn-89  14:4  4 


load 


.PBDBIEM:  'SIAT' 

H/P  SENSITIVITY: 
-208.G2dB  (vclts/uPa) 

.DISPLACEMENT  SCALE  R 

i  u.ijyu2byt;-uj  m — I 


.DISPLACEMENT  SCALE  Z 

I  0.13902SgE-03  m  I 

RMIN=  -0.0020000  m 
RMAX=  0.1020000  m 
ZMIN=  -0.0520000  m 
ZMAX”  0.0520000  m 
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8  SERIES  5;  ANALYTICAL  FLUID  PROBLEM 

A  cross  sectional  diagram  of  the  structure  modelled  in  this  series  is  ^ven  in 
Figure  8,0.1.  It  is  a  10  cm  diameter  piston  that  vibrates  axially  in  an  infinite 
hard  baffle.  The  half  space  over  the  piston  and  baffle  is  sea  water.  The  outer 
radius  of  the  sea  water  sphere  used  in  the  analysis  is  0.165  m. 

The  mesh  of  the  F.E.  model  with  node  numbers  is  shown  in  Figure  8.0.2,  while 
Figure  8.0.3  gives  the  element  numbering.  The  piston  is  modelled  with  six 
SHELL  (Type  10)  elements. 

The  model  shown  in  Figures  8.0.2  and  8.0.3  has  been  changed  from  that  in  the 
original  Examples  Manual  [1.1c].  The  main  change  has  been  to  reduce  the 
number  and  size  of  the  fluid  elements  to  increase  accuracy  and  to  reduce  the 
rim  time.  Testing  that  resulted  in  these  changes  is  reported  in  Ref.  [4.2] 


Figure  8.0.1:  Series  5  Problem  -  Diagram  of  Piston  in  Infinite  Baffie 

az 
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8.1  PROBLEM  15 

Analysis  Type:  DRIVE  -  piston  in  infinite  hard  baffle 

Input  File:  D1 5  .DAT 

The  piston  is  driven  at  a  vibration  amplitude  of  1  mm  by  appl3dng  fixities  in  the 
Z  direction  on  all  H  nodes  of  the  SHELL  elements.  Only  one  excitation 
frequency  of  12,000  Hz  is  used  in  the  analysis. 

Test  results  reported  in  [4.2]  show  that  this  model  is  acoustically  accurate  at 
frequencies  up  to  18,000  Hz.  It  was  shown  that  FLUID  element  dimensions 
should  not  exceed  0.4  wavelength  at  the  highest  frequency  where  accurate 
results  are  sought. 

A  plot  of  far-field  directivity  at  12,000  Hz  is  shown  in  Figure  8.1.1.  The  peak 
pressure  of  269.28  dB  reljjPa  @lm  is  the  same  as  predicted  by  anal3dical  theory 
[4.2].  When  an  acoustic  radiation  problem  is  analysed  as  a  hemisphere,  a 
mirror  of  symmetry  is  assumed,  hence  the  directivity  pattern  is  repeated  for  the 
lower  hemisphere.  Under  these  conditions,  a  hard  baffle  condition  on  the 
symmetry  plane  is  automatic.  A  soft  baffle  could  be  created  by  imposing 
pressure  fixities  of  zero  on  all  fluid  nodes  on  the  plane,  but  this  soft  baffle  would 
extend  only  to  the  spherica’.  boundary  of  the  near-field  region  analysed. 

Figure  8.1.2  shows  near-field  pressure  contours  at  12,000  Hz. 
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Figure  8.1.1:  Problem  15  -  Piston  Directivify  at  12,000  Hz 


PEOBLEK  #15  -  Piston  in  Hard  Baffle 
269.30 


DIRECTIONAL  RESPONSE 
(dB  re  1  ^iPa/volt  8  1  m) 

•DRIVE  FREOJENCY: 

1200  0.0  Hz 

-DIRECTIVITY  INDEX: 

-11.07  64  dB 

-PEAK  TRANSMIT  RESPONSE: 
269.277  dB 

.RECEIVING  RESPONSE  0  EEG; 

0.00 


ll-DEC-89  12:54 


Figure  8.1.2:  Problem  15  -  Near-Field  Pressure  Contours 

at  12,000  Hz 


PROBLEM  #15  -  Piston  in  Hard  Baffle 
ENTIRE  REGION  RADIUS-0.1650 
ISCAMPLITUDE  OONTOUPS 


FREQUENCY:  1200  0.  00  Hz 

_  _  POSITIVE 

—  ZERO 

.  NEGATIVE 

ISOAMPLITUDE  CCNTOURS 
.MAXIMJM  AMPIIIUDE: 

0.2286280E+09  Pa 
.MINIMJM  AMPLITUDE: 
0.2084398E+08  Pa 


.REFERENCE  AKPLITCDE: 

0  .  690  3  2  7  2E+08  Pa 
.AMPLITUDE  CCNTOURS: 


10.40000 

dB 

9.100000 

dB 

7.800001 

dB 

6.500000 

dB 

5.200000 

dB 

3.90  0  00  0 

dB 

2 . 60  0  0  0  0 

dB 

1.300000 

dB 

O.OOOOOOQE+00  dB 

-1.30  0  00  0 

dB 

-2.60  0  0  0  0 

dB 

-3.900000 

dB 

-5.200000 

dB 

-6.500000 

dB 

-7.800001 

dB 

-9.100000 

dB 

-10.40000 

dB 

ll-DEC-89 

13:07 
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9  SERIES  6:  COMPLIANT  FLUID 

This  series  demonstrates  the  use  of  Compliant  Fluid  modelled  with  SOLID 
elements  to  simulate  a  20-cm  diameter  spherical  air  bubble  pulsating  in  sea 
water.  Figure  9.0.1  is  a  cross  sectional  diagram  of  the  structure  analysed. 

The  mesh  of  the  F.E.  model  is  shown  in  Figure  9.0.2.  The  material  properties 
of  the  air  within  the  bubble  are  determined  assxnning  that  the  bubble  is  placed 
at  10  metres  depth,  where  the  absolute  hydrostatic  pressure  in  the  water  is  two 
atmospheres.  Values  for  density  and  bulk  modulus  are  2.585  kg/m3  and 
0.285E+6  Pa,  respectively.  The  adiabatic  bulk  modulus  of  the  air  is  used,  as  this 
is  a  d3mamic  problem;  for  static  problems,  the  isothermal  bulk  modulus  would 
be  used. 

Figures  9.0.2  (a,b,c)  and  9.0.3  (a,b)  give  details  of  node  and  element  numbering. 
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Figure  9.0.1:  Series  6  Problems  -  Air  Bubble  in  Water 


Figure  9.0.2  (a):  Series  6  Problems  -  Outer  Fluid  Node  Numbering 

Series  6,  Spherical  air  bubble 


n  o 
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Figure  9.0.3  (a);  Series  6  Problems  -  Outer  Element  Numbering 


Series  S,  Spherical  air  bubble 
GEOyETRy  PIOTTINS 


HjEMENT  NIMBERS 

.D3MENSICNS: 
(Hsnrra:  resecm 
PMIN=  -O.^QTOO 
PiyK«=  1.0170001 
ZM[N=  -O.Q17QQQ0 
ziyw#=  1.0170001 

.MKTERIRL  TIPES: 

1 

■ETFJ^EMT  TYIES: 
4-7-15 


y  srv/  4o\ 

28' 


I'O  15-EEC-89  15:31 


Figure  9.0.3  (b):  Series  6  Problems  -  Inner  Element  Numbering 

Series  6,  Spherical  air  bubble 


««««  o  sees 
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9.1  PROBLEM  16 

Analysis  Type:  DRIVE  -  spherical  air  bubble 

Input  File:  D16.DAT 

Since  only  the  upper  hemisphere  of  the  complete  structure  is  analyzed,  all  "air" 
nodes  on  the  symmetry  plane  must  be  restrained  in  the  Z  direction. 

A  imiform  volumetric  harmonic  body  force  was  applied  to  the  SOLID  elements 
comprising  the  air  by  placing  thermal  loads  on  the  elements. 

Frequency  response  plots  of  radiated  far-field  pressure  (dB  re  IpPa  @  Im)  and 
radial  displacement  of  the  surface  of  the  bubble  at  Node  69  are  given  in  Figure 
9.1.1.  The  resonance  frequency  is  45.934  Hz  with  a  Q  of  51.57.  The  resonance 
frequency  predicted  by  analytical  theory  [4.2]  is  45.9  Hz  with  a  Q  of  52.0. 

The  pressure  in  the  air  bubble,  given  as  the  pressure  in  the  stress  printout  for 
each  element,  is  constant  at  11000  Pa.  This  is  equal  to  the  nodal  pressure  in 
the  water  at  the  surface  of  the  bubble. 
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10  SERIES  7:  RING 

To  demonstrate  the  use  of  RING  elements,  the  same  tubular  shell  structure  is 
used  as  in  Problem  13,  with  RING  stiffeners  added.  A  representative  cross 
sectional  diagram  is  given  in  Figure  10.0.1.  The  node  munbering  is  the  same 
as  in  Figure  7.0.2.  Ten  RING  elements  (#11  to  #20)  are  added  to  the  structure  of 
Figure  7.0.3,  and  are  located  at  the  odd  numbered  H  nodes  of  the  shell.  The 
cross  section  of  a  single  ring  is  5  mm  high  and  2  mm  thick.  The  end  nodes,  1 
and  21  are  stiffened  with  single  ring  ribs,  while  the  rest  of  the  appropriate 
nodes  receive  double  reinforcing. 


Figure  10.0.1:  Series  7  Problems  -  SHELL  Tube  with  RING  Stiffening 


(cm) 
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10.1  PROBLEM  17 

Analysis  Type:  EIGEN  -  RING-stififened  SHELL  tube, 

(a)  with  no  additional  inertia, 

(b)  with  nodal  masses. 

Input  Files:  D17A.DAT 

D17B.DAT 

In  Case  (a),  the  original  structure  of  Figure  7.0.1  (Problem  13)  is  stiffened  with 
alumintun  RING  elements  added  to  odd  numbered  H  nodes  of  the  shell.  No 
inertia  term  (density)  is  specified  in  the  material  properties  data  for  the  ring,  so 
that  the  stiffness  alone  may  be  modelled.  The  end  nodes,  1  and  21  are  free  to 
move  only  in  the  R  direction,  corresponding  to  the  condition  of  Figure  7.1.2. 

The  parameter  POLANG  in  the  element  data  for  RING  elements  defines  the 
number  of  elementary  rings  assigned  to  any  given  RING  element;  thus, 
POLANG  is  unity  for  elements  1  and  10,  and  2  for  the  remaining  RING 
elements. 

In  Case  (b),  nodal  masses  corresponding  to  the  masses  of  each  respective  RING 
element  are  specified  at  the  appropriate  nodes.  This  essentially  takes  the  place 
of  the  density  term  that  could  have  been  included  in  the  material  data  for  the 
RING,  Note  that  the  total  mass  of  the  ring  is  added  to  each  DOF  on  which  it  is 
effectively  acting.  If  the  inertial  moments  of  the  RING  elements  are  to  be 
modelled,  the  density  must  be  included  in  the  material  data,  as  the  added 
masses  cannot  represent  inertial  moments. 

In  both  cases  the  radial  displacement  of  Node  1  is  used  as  the  normalizing  DOF 
and  is  set  to  1.  The  values  of  ERVAL  and  ERVEC  are  set  to  1.0E3— 8  and  l.OE— 6, 
respectively. 

The  results  show  the  first  natural  frequency  of  the  RING- stiffened  shell  to  be 
21,756  Hz  without  the  nodal  masses  (Case  (a)).  After  adding  nodal  masses,  the 
frequency  drops  to  15,418  Hz  (Case  (b)).  The  mode  shapes  are  both  essentially 
the  same,  and  are  the  same  as  shown  in  Figure  7.1.2  for  Problem  13.  As  in 
Problem  13,  the  pure  hoop  mode  is  again  Mode  4,  but  its  fi*equency  is  lowered  to 
16,933  Hz. 
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11  SERIES  8:  DAMPING  APPLICATIONS 

The  model  used  in  Problem  9  (Figure  5.0.1)  was  chosen  for  demonstration  of 
damping  employment  in  solids  and  fluids.  Node  numbering  is  the  same  as  in 
Figure  5.0.3  (a,b).  Element  numbering  for  Problem  18A  is  the  same  as  in 
Figure  5.0.2.  For  Problem  18B,  an  additional  five  VISCous  fluid  elements  (#81 
to  #85)  are  added  to  the  structure,  using  existing  fluid  nodes. 


11.1  PROBLEM  18 

Analysis  Type:  DRIVE  -  damping  in  the  PT  solid  ring  (a), 

-  damping  in  a  VISCous  fluid  layer  (b) 

Input  Files:  D1 8A(1 ,2).DAT 

D18B(1,2,3,4).DAT 

In  Case  (a)  the  stiffness  and  mass  distributed  damping  coefficients,  |j.k  and  p-m 
equal  to  0.03  (3%),  were  arbitrarily  chosen.  Otherwise  the  data  did  not  differ 
from  those  used  in  Problem  9.  The  two  damping  coefficients  were  applied 
separately,  and  resonance  searches  were  applied  in  both  cases  in  the  range 
from  250  to  450  Hz.  Transmitting  voltage  responses  for  the  two  types  of 
damping  are  plotted  in  Figure  11.1.1,  together  with  the  result  of  Problem  9  for 
no  damping.  Corresponding  efficiency  plots  are  given  in  Figure  11.1.2  for  the 
two  damped  cases.  (The  efficiency  for  the  imdamped  structure  is  100  percent.) 

In  Case  (b),  the  material  damping  coefficients  for  the  PT  ring  were  set  to  zero, 
and  a  viscous  dissipating  fluid  layer  was  inserted  on  the  ring.  This  comprises 
five  Vise  elements,  laid  on  F  nodes  of  the  FTOS  elements. 

The  appropriate  material  properties  matrix  for  the  VISC  elements  was  added 
to  the  data.  Viscous  damping  is  invoked  by  setting  the  "density"  to  1.0  for 
pressure  var3dng  damping  and/or  setting  the  stiffness  coefficient  cn  to  1.0  for 
pressure  gradient  varying  damping;  these  are  flags  and  have  no  relevance  to 
their  normal  properties.  The  magnitude  of  the  damping  is  applied  in  the  same 
way  as  for  solid  material  damping:  |ik  for  stiffness  distributed  damping  and  |j.m 
for  mass  distributed  damping,  although  these  parameters  are  no  longer 
fractions  between  zero  and  one.  The  four  types  of  viscous  damping  were 
applied  separately  to  the  model  and  a  resonance  search  was  done  in  the  same 
range  as  for  material  damping.  Values  of  |j,k  and  |j.ni  were  selected  so  as  to 
yield  efficiencies  near  80%  at  the  375  Hz  resonance. 

Plots  of  the  transmitting  voltage  responses  with  viscous  damping  are  compared 
to  those  for  the  undamped  structure  in  Figures  11.1.3  to  11.1.6.  The 
corresponding  efficiencies  are  plotted  in  Figure  11.1.7. 
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Figure  11.1.1:  Problem  18(a)  -  Transmitting  Voltage  Responses 

of  a  Free-Flooding  Ring  with  DLEferent  Material 
Damping 


PKDBLEM  #18A,  PT  Piezoelectric  Ring  with  Material  Damping 
— a —  Stiffness  damping  .  Mass  damping  No  damping 


FREQDENCY  <Hz)  2-JAN-90  15:40 


Figure  11.1.2:  Problem  18(a)  -  Efficiencies  with  Material  Damping 


PROBLEM  #18A,  PT  Piezoelectric  Ring  with  Etiterial  Damping 

-  Stiffness  damping,  3%  • — Mass  damping,  3% 
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Figure  11.1.3:  Problem  ISOb)  -  Transmitting  Voltage  Response  for  a 

Free-Flooding  Ring  with  Viscous  Pressure  Damping, 
|ik  =  2.0E+8 


PROBLEM  #18B,  PT  Piezoelectric  Ring,  VISC  Fluid  layer.  Press .K 
— a —  Press.  K  damping  — K--  No  damping 


Figure  11.1.4:  Problem  18(b)  -  Transmitting  Voltage  Response  for  a 

Free-Flooding  Bing  with  Viscous  Pressure  Damping, 
!Xm=  30 


PRCBIEM  #18B,  PT  Piezoelectric  Ring,  VISC  Fluid  layer.  Press. M 
-  Press .  M  damping  -  No  damping 
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Figure  11.1.5:  Problem  18(b)  -  Trajosmittiiig  Voltage  Response  for  a 

Free-Flooding  Ring  with  Pressme  Gradient  Damping, 
{i-k  =  1.5E+7 


PROBLEM  #18B,  PT  Piezoelectric  Ring,  VISC  Fluid  Layer,  P.grad.K 
— S —  Pr.grad.  K  damping  — K--  No  damping 


Figure  11.1.6:  Problem  18(b)  -  Transmitting  Voltage  Response  for  a 

Free-Flooding  Ring  with  Pressiire  Gradient  Damping, 
jim  *  3 


PROBLEM  #18B,  PT  Piezoelectric  Ring,  VISC  Fluid  Layer,  P.grad.M 
-  Pr.grad.  M  damping  -  No  damping 
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Figure  11.1.7:  Problem  18(b)  -  Efficiencies  for  a  Free-Flooding  Bing 

with  Different  l^pes  of  ^^cous  Fluid  Damping 

PBOBIEM  #18B,  PT  Piezoelectric  Ring,  VISCous  fluid  layer 
-  Pr.  K  .  Pr.  M  - P.Grad  K  .  P.Grad  M 


We  note  that  the  solid  material  stiffness  damping  has  a  profound  effect  on  the 
resonance  frequency  of  a  free  flooding  ring,  while  the  other  forms  of  damping 
have  virtually  no  effect  on  the  frequency.  In  all  cases,  the  dependence  of 
efficiency  on  frequency  is  greater  for  stiffness  damping  than  for  mass 
damping.  Further  investigation  of  damping  in  a  free  flooding  ring  transducer 
will  be  reported  in  a  future  DREA  Note. 
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12  SERIES  9:  RING  ARRAY,  RIGID  NODE 

Two  PT  rings,  each  made  of  48  staves  of  Channel  5400  ceramic,  are  immersed 
in  sea  water  and  separated  by  an  aluminum  ring  as  shown  in  Figure  12.0.1. 
The  mesh  of  the  F.E.  model  is  given  in  Figure  12.0.2.  Figures  12.0.3  (a-d)  and 
12.0.4  (a,b)  give  details  of  node  and  element  munbering,  with  the  radial , 
dimensions  expanded  relative  to  the  axial  dimensions.  The  E  type  Nodes  975 
and  987  coincide  with  F  type  Node  376,  and  are  placed  at  the  origin  of  the  system 
of  coordinates. 

•  For  this  edition  of  the  Examples  Manual,  the  model  was  changed  by  deleting 
the  four  outer  layers  of  fluid  elements  from  the  original  model.  This  reduces 
the  size  of  the  problem  with  a  possible  improvement  in  accuracy  as  noted  under 
Problem  15.  Also,  the  eight  elements  representing  the  two  PT  rings  have  been 
placed  in  two  element  groups,  one  for  each  ring  (NELG  =  2  on  Card  10).  When 
all  eight  elements  are  plac3d  in  one  group,  the  radiated  power,  admittance  and 
efficiency  are  in  error. 
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Figure  12.0.1:  Series  9  Problems  -  Cross-Sectional  Diagram  of  Model 

^2  (cm) 
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Figure  12.0.2:  Series  9  Problems  -  F.E.  Mesh 


CCMPIETE  GRID 

. DIMENS ICNS 
RMI»=  -0.1024000M 
RMAX=  0.3024001M 
ZMIN=  -0.2024000M 
ZMAX-  0.2024000M 


.MATERIAL  TYPES: 
1-2-4 


■ELEMENT  TYPES: 
4-6-7-12-14-15 


22-JAN-90  13:22:26 


Figure  12.0.3  (a):  Series  9  Problems  -  Outer  Node  Numbering 


NDDE  NUMBERS 

■  DIMENSICNS: 
(ENTIRE  REGION) 
RMIN=  -0.0100000 
RMAX=  0.2100000 
ZMIN=  -0.2027520 
ZMAX=  0.2112000 

■MATERIAL  TYPES: 

1 


■  ELEMENT  TYPES: 
4-6-7-12-14-15 


5-FEB-90  13:38 


Series  9,  Two  Fr^  PT  Rings  with  SOLID  Spacer 
GEOMETRY  PLOTTING 


Series  9,  Two  Free  PT  Rings  with  SOLID  Spacer 
GECMETRY  PLOTTING 
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Figure  12.0.3  (b):  Series  9  Problems  -  Middle  Node  Numbering 


Series  9,  Two  Free  PT  Rings  with  SOLID  Spacer 
S:OtETRY  PLOTTIMS 


NODE  NUMBERS 

•  DIMENSIONS: 

RMIN=>  -0.0023760  m 
RMAX=  0.1163  7  65  m 
ZMIN”  -0.10  97  627  ra 
Zmx-  0.1137714  m 

.MATERIAL  TYPES: 

1 

.EIEMENT  TYPES: 
4-6-7-12-14-15 


5-FEB-90  14:03 


Figure  12.0.3  (c):  Series  9  Problems  -  Inner  Fluid  Node  Numbering 


Series  9,  Two  Free  PT  Rings  with  SOLID  Spacer 


NODE  NUMBERS 

DIMENSIONS: 

RMIN=  0.0332498  m 
RMAX=  0.0688755  m 
ZMIN=  -0.0561145  m 
ZMAX=  0.0578879  m 

.MATERIAL  TYPES: 

1 

.ELEMENT  TYPES: 
4-6-7-12-14-15 


5-FEB-90  14:05 
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Figure  12.0.3  (d):  Series  9  Problems  -  Solid  Node  Numbering 


Series  9,  Two  Free  PI  Rings  with  SOLID  Spacer 


GEOMETRY  PLOTTING 


NODE  NIMBERS 

.DIMENSIONS 
RMIN=  0.0235839M 
RMAX=  0.0801373M 
2MIN=  -0.0607200M 
ZMAX=  0.0607200M 

•  MATERIAL  TYPES: 
2-4 

. ELEMENT  TYJES : 
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Figure  12.0.4  (a):  Series  9  Problems  -  Outer  Element  Numbering 
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Figure  12.0.4  (b):  Series  9  Problems  -  Inner  Element  Numbering 
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12.1  PROBLEM  19 

Analysis  Type:  DRIVE  -  two  free  flooding  PT  rings  with  SOLID  spacer 

Input  Files:  D1 9  .DAT 

The  model  shown  in  Figures  12.0.3  and  12.0.4  is  analyzed.  A  6000  Hz,  1-V 
amplitude  excitation  is  applied  to  both  E  type  nodes.  The  inner  midplane  nodes 
of  each  PT  ring  and  the  spacer  are  restrained  in  the  Z  direction,  i.e..  Nodes  185, 
458,  and  324. 

Plots  of  displacement,  far-held  directivity,  and  near-held  pressure  contours  are 
shown  in  Figures  12.1.1  through  12.1.3,  respectively. 


Figure  12.1.1:  Problem  19  -  Displacement  at  6000  Hz  with  SOLID  Spacer 


PRDBI£M  #19  Two  Free  PT  Rings  vdth  SOLID  Spacer 
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Figure  12.L2:  Problem  19  -  Fai>Field  Directivify  at  6000  Hz 

with  SOLID  Spacer 
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Figure  12.1.3:  Problem  19  -  Near-Field  Pressure  Contours 

with  SOLID  Spacer 
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12.2  PROBLEM  20 

Analysis  Type:  DRIVE  -  two  free  flooding  PT  rings  with  a  RIGID  node 

replacing  the  spacer 

Input  FHe:  D20.DAT 

Using  the  model  from  Problem  19,  the  aluminum  spacer  ring  is  replaced  by  a 
single  R  type  (rigid)  node.  Of  all  the  S  nodes  of  the  spacer,  only  the  central  Node 
282  is  retained  and  changed  to  R  type.  The  R  node  becomes  a  common  solid 
node  for  all  the  FTOS  elements  surrounding  the  spacer,  while  the  F  nodes 
remain  the  same  as  in  Problem  19.  Again,  the  same  nodes  of  the  PT  rings  as 
well  as  the  R  node  are  restrained  in  the  Z  direction.  The  same  excitation  as  in 
Problem  19  is  used. 

Results  of  the  analysis  are  practically  the  same  as  in  Problem  19.  Figure  12.2.1 
shows  a  plot  of  near  field  pressure  contours  for  comparison  to  Figure  12.1.3. 


Figure  12.2.1:  Problem  20  -  Pressure  Contoiu^  with  a  Rigid  Node 

as  the  Spacer  Ring 
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13  SERIES  10:  TORSIONAL  PROBLEMS 

An  aluminum  cylinder  20  cm  in  diameter  and  40  cm  long  is  used  to 
demonstrate  torsional  problems.  A  diagram  of  the  structure  is  given  in  Figure 
13.0.1.  The  Finite  Element  mesh  is  shown  in  Figure  13.0.2  with  element 
numbers,  and  node  numbers  are  given  in  Figure  13.0.3. 

Since  real  torsional  vibration  modes  cannot  couple  to  a  fluid  acoustically,  fluids 
will  never  be  used  meaningfully  in  a  torsional  problem. 

Driving  of  torsional  deformations  piezoelectrically  has  not  yet  been 
demonstrated  with  MAVART.  This  will  require  axial/radial  poling  with  a 
tangential  driving  field  or  tangential  poling  with  axial/radial  driving  fields, 
implying  that  some  modification  to  PAR  and/or  PT  elements  may  be  necessary. 


Figure  13.0.1:  Series  10  Problems  -  Solid  Cylinder  Model 

for  Torsional  Problems 


(cm) 
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Figure  13.0.2:  Series  10  Problems  -  Element  Numbering 


Series  10,  Torsional  m=0  Examples 
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Figure  13.0.3:  Series  10  Problems  -  Node  Numbering 
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13.1  PROBLEM  21 

Analysis  Type:  STATC  -  Torsional  load 

Input  FHe:  D21.DAT 

Torsional  analysis  is  invoked  by  setting  MSYM  =  -1  in  Card  2,  with  the  Fourier 
mode  number  IMODE  =  0.  The  nodes  on  one  end  of  the  cylinder  are  fixed  in  R,  Z, 
and  0,  and  a  100  N  torsional  load  is  applied  at  Node  25  on  the  other  end  of  the 
cylinder.  As  this  load  is  at  a  radius  of  0.1  m,  the  torque  is  10  Nm. 

From  simple  elastic  theory,  the  displacement  in  radians  for  a  circular  cylinder 
of  radius  a  and  length  L,  subjected  to  a  torque  T  is  given  by 

0  =  2TL  I  (.Gna^)  (13.1) 

where  G  =  0.2664E11  Pa  is  the  shear  modulus  of  the  aluminum.  Eq.  13.1 
predicts  a  tangential  displacement  of  0.9559E— 6  radian  at  the  end  of  the 
cylinder  for  a  uniform  torque  of  10  Nm.  MAVART  predicts  a  displacement  of 
0.9927E-6  rad  at  Node  25,  and  0.9152E-6  rad  at  Node  23,  which  is  at  midradius. 
This  variation  is  due  to  the  load  being  applied  on  the  outer  radius,  and  not 
uniformly  across  the  end  of  the  cylinder. 

The  graphics  program  GRAFl  does  not  support  plotting  of  torsional 
displacements,  hence  no  graphic  output  is  displayed. 


13.2  PROBLEM  22 

Analysis  Type:  EIGEN  -  Torsional  modes 

Input  File:  D22.DAT 

A  modal  analysis  was  conducted  on  the  model  depicted  in  Figure  13.0.1.  The 
same  zero  fixities  were  applied  and  the  100  N  load  was  removed,  so  the 
vibrating  structure  is  a  fixed-free  cylinder.  Torsional  modes  will  occur  when 
the  cylinder  length  is  an  odd  multiple  of  one  quarter  wavelength.  The  shear 
velocity  is  given  by  Cg  =  V(G/p ),  where  p  =  2700  kg/m^  is  the  density  of  the 
aluminum.  The  calculated  fundamental  frequency  for  the  0.4  m  long  cylinder 
is  1963.1  Hz.  MAVART  predicts  1963.7  Hz,  in  very  close  agreement. 
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14  ST^RTPIS  U  PROBLEMS  -  NONLINEAR  EFFECTS 

A  cylinder  as  shown  in  Figure  14.0.1  is  used  to  demonstrate  effect  of 
geometrical  nonlinearity  on  longitudinal  vibrations  of  a  umform  bar. 

The  finite  element  mesh  consists  of  16  quadrilateral  solid  elements  (Type  12) 
and  85  nodes. 

Figure  14.0.1  Series  11  Problems  -  Solid  Cylinder  Model  for 

Nonlinear  Effects 

I  P 


10  cm 
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The  free  end  of  the  circular  cylinder  is  subjected  to  uniform  external 
compression  of  magnitude  p.  The  other  end  of  the  cyhnder  is  fixed  in  R,  Z  and  (j). 


14.1  Problem  23 

Analysis  Type:  Static-nonlinear 

Input  File:  d23.dat 


A  nonlinear  effect  was  analyzed  applying  several  pressure  values  to  the  free 
end  of  the  cylinder.  Resxilting  displacement  of  the  cylinder  end  at  the  centre 
(R=0)  is  shown  in  Table  14.01. 


Table  14.01 

Cylinder 

Displacement 

p  (Pa) 

Linear 

Deformation 

Nonlinear 

Deformation 

Compression/Tension 

Linear 

Theory 

106 

-5.4830.10-6 

-5.4863.10-6/-5.4797.10-6 

-5.5811.10-6 

107 

-5.4830.10-5 

-5.51 60.1 0-5/-5.4505.1 0-5 

-5.5811.10-5 

108 

-5.4830.10-4 

-5.8379.1 0-4/-5.1 783.1 0-4 

-5.5811.10-4 

14.2 

Problem  24 

Analysis  Type:  EIGEN  -  Nonlinear  effects 

Input  File:  d24.dat 


A  modal  analysis  was  conducted  on  the  model  shown  in  Figure  14.01.  The 
same  zero  fixities  were  applied  and  the  pressure  load  was  removed.  Additional 
stiffness  matrices  resulting  from  the  pressure  load  were  generated  by  running 
the  STATIC  problem  with  the  input  NLGEO  flag  set  to  1 ,  prior  to  running  a 
corresponding  EIGEN  problem.  For  the  additional  stiffness  to  be  included  in 
the  EIGEn  analysis,  the  input  NGLEO  flag  was  set  to  2.  Results  of  the 
nonlinear  EIGEN  analysis  are  summarized  in  Table  14.02. 
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Table  14.02  Effect  of  additional  stiffness  on  Natural  Frequency  of  the 


p 

Linear 

Problem 

Cylinder  (Hz) 

Nonlinear 

Problem 

C  ompression/Tension 

Theory 

0 

3266.4 

3266.4 

3220.1 

106 

3266.4 

3263.2/3269.6 

107 

3266.4 

3234.6/3297.9 

108 

3266.4 

2929.7/3565.9 

The  wave  longitudinal  velocity  is  C^s/Wp  where  p=2700  kg/m^  is  the  density  of 
almninium  and  E=7.167.10l0  is  the  Young's  modulus.  The  calculated 
fundamental  frequency  of  a  cylinder  of  length  L=0.4  m,  is  f  =  Co/4L  =  3220.1  Hz. 
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15  SUMMARY 

1.  With  few  exceptions,  all  example  problems  from  the  Examples  Manual 
[l.lc]  have  been  run  successfully  at  DREA,  using  one  or  the  other  of  the  latest 
two  versions  of  MAVART  (MAVART  11.0  and  12.0)  .  Some  of  the  problems' 
input  data  required  modification  or  minor  correction  to  3deld  satisfactory 
results. 

2.  The  FTOS  Reversibility  deficiency  in  MAVART  has  been  resolved,  and 
Problem  16  has  been  reformulated  with  Compliant  Fluid  modelling  in  keeping 
with  the  findings  of  [4.4]. 

3.  A  deficiency  in  MAVART  remains: 

•  It  is  not  clear  whether  piezoelectric  drive  of  torsional  deformations  can  he 
acommodated  without  some  modification  to  the  PAR  and/or  PT  elements 
(Series  10  problems). 

4.  The  postprocessor  program  GRAFl  has  been  used  to  display  a  broad 
range  of  analysis  data  from  all  example  problems  except  the  torsional 
problems.  As  well,  GRAFl  has  been  used  to  display  model  geometry  for  all 
examples. 
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16  RECOMMENDATIONS 

1.  An  example  problem  that  demonstrates  the  piezoelectric  drive  of 
torsional  deformations  should  be  added  to  the  Series  10  Problems.  If  necessary, 
MAVART  should  be  modified  to  allow  such  a  driving  mode. 

2.  The  DREA  postprocessor  GRAFl  should  be  updated  to  support  plotting  of 
tangential  displacements  for  both  torsional  problems  and  for  problems  where 
the  Fourier  symmetry  m  is  greater  than  zero. 


16-1 


P153521.PDF  [Page:  118  of  124] 


P153521.PDF  [Page:  119  of  124] 


17  REFEDEIENCES 

[1.1]  (a)  User's  Manual  for  Program  MAVART, 

(b)  Theoretical  Manual  for  Program  MAVART, 

(c)  Examples  Manual  for  Program  MAVART  Vol.  I  &  11, 

DREA  Contractor  Report  CR/87/442,  E.  L.  Skiba,  June,  1987. 

[1.2]  User's  Manual  for  Program  MAVART,  E.  L.  Skiba  and  G.  W.  McMahon, 
DREA  Technical  Memorandum  TM  91/201,  January,  1991. 

[1.3]  Examples  Manual  for  Program  MAVART,  G.  W.  McMahon  and 

E.  L.  Skiba,  DREA  Technical  Memorandum  TM  90/204,  August,  1990. 

[4.1]  Piezoelectric  Technology,  Data  for  Designers,  Clevite  Corporation, 
Piezoelectric  Division  (1965). 

[4.2]  MAVART  Testing  -  Part  I,  G.  W.  McMahon,  DREA  Note  SP/89/8, 

July  1989. 

[4.3]  Research  and  Development  of  Spherical  Transducers,  Final 
Development  Report,  H.  Jaffe,  F.  Rosenthal,  and  H.  Baerwald,  Clevite 
Research  Center,  ONR  Contract  Nonr  1825(00),  April  1957. 

[4.3]  Resolution  of  FTOS  Reversibility  Issue,  DREA  Contractor  Report 
CR/91/467,  E.L.  Skiba,  June  1991. 


17-1 


P153521.PDF  [Page:  120  of  124] 


P153521.PDF  [Page:  121  of  124] 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 
[highest  dasslflcallon  of  Title,  Abstract,  Keywords] 


DOCUMENT  CONTROL  DATA 

(bounty  ctassfficaticm  of  tiUd,  body  of  abstr^  and  indexing  annotation  must  be  ent^ed  wh^  tiie  overall  dooim^t  is  classified) 

1.  UHttiiNAiUH  (the  name  and  adoressot  the  organization  prepanng  me  oocumera 

Organizations  for  whom  the  document  was  prepared,  e.g.  Establishment  sponsoring 
a  tantractoris  report,  or  tasking  agency,  are  entered  in  section  8.) 

Acres  International  Limited 

5259  Dorchester  Rd.,  Niagara  Falls, 

Ontario,  Canada,  L2E  6W 1 

2.  SECUnn  y  CLASSIHCATtON - 

(overall  seojifty  ciassificaUon  of  the  document 
including  special  warning  tetms  If  appCcable). 

Unclassified 

3.  n  1  Lh  (me  (^mpiete  document  ti&e  as  indicated  on  me  title  page,  its  dassitication  snould  be  indicated  by  me  appropnate 

abbreviation  (S,C,R  or  U)  in  parentheses  after  the  title). 

Examples  Manual  for  Program  MAVART 

4,  AUfHOFis  (Last  name,  tirst  name,  middle  initial.  It  military,  show  rank,  e.g.  uoe,  Maj.  John  b.) 

E.  L.  Skiba,  G.  W.  McMahon  and  Z.  Wozniak 

DATE  OF  PUBUCATION  (month  and  year  of  publication  ot 
document) 

January  1995 

6a.  NO  OF  PAGES  (total 
containing  infoimatlon  Include 
Annexes,  Appendices,  etc). 

91 

6b.  NO.  OF  REFS  (total  cited  in 
document) 

7 

/.  utSCRlP  1  ivt  NO  1  ta  (me  category  ot  me  document,  e.g.  technical  report,  technical  note  or  memorandum,  it  appropr 
progress,  summary,  armual  or  final.  Give  the  inclusive  dales  when  a  specific  reporting  period  is  awered). 

Contractor  Report 

ate,  enter  the  type  ot  tepoil,  e.g.  intenm. 

8.  SPUNaUHiNU  AG  1 1 VI 1 Y  (the  name  ot  the  department  project  ottice  or  laboratory  sponsonng  me  research  and  devmopment.  include  me  address). 

Defence  Research  Establishment  Atlantic,  9  Grove  St.,  PO  Box  1012, 

Dartmouth,  NS,  Canada,  B2Y  3Z7 

»a.  RHUJtu  1  UM  UHAN 1  NU.  (11  appiopnate,  me  appiicaoie  researcn 
and  developmojt  project  or  grant  number  under  which  the  document 
was  written.  Please  specify  whether  prefect  or  grant). 

lAU-13 

yo.  UON 1  HAG  1  NU.  (It  appropnate,  me  applicable  number  under 
which  the  document  was  written). 

W7707-4-2206/0 l-OSC 

10a.UHIGINA10H'tS  DOCUMEN 1  NUMBER  (the  ollidal  document 
number  to  which  the  document  is  identified  by  the  originating 
activity.  This  nunrtjer  must  be  unique  to  this  document). 

DREA/CR95/416 

TT!  UUUUMtN  1  AVAlLABILl  I Y  (any  limitations  on  further  dissemination  of  the  document,  other  than  those  imposed  oy  security  classilication) 


(  X )  Unlimited  distribution 

I  )  Distnbutiwi  limited  to  defence  departments  and  defence  contractors;  further  distnbuBon  only  as  approved 
{  )  Distritwtion  limited  to  defence  departments  and  Canadian  defence  contractors;  further  distribution  only  as  approved 
(  i  Distribution  limited  to  govemmait  departments  and  agencies;  further  distribution  only  as  approved 
(  j  Dislribulion  limited  to  defence  departments;  further  distribution  only  as  approved 
(  )  Oiher  (please  specify): 

TS:  DOCUMfcN  I  ANNUUNCtrvitN  I  (any  limitation  b  the  bibliographic  announcement  ot  this  document.  I  nis  will  normally  conespond  to  the  Document  Availability 
(11).  However,  where  further  distribution  (beyond  the  audience  spealied  in  11)  is  possible,  a  wider  announcement  audience  may  be  selected). 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 

DCD03  2/06m7-M 


P153521.PDF  [Page:  122  of  124] 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 

ABB !  HAU I  la  cinet  and  taouial  summaiy  oi  tne  document,  it  may  also  appear  elsewhere  in  me  Body  ot  me  document  itsell.  It  is  highly  desirawe  that  the  abstract 
of  classified  docutaenls  be  unclassified.  Each  paragraph  of  the  abstract  shall  begin  with  an  indication  of  fhe  security  dassification  of  the  infonmation  in  the  paragraph 
{unless  the  document  Itself  is  undassified)  represented  as  (S),  (C),  (R),  or  (U).  It  is  nra  necessary  to  include  here  abstracts  in  both  official  languages  unless  the  text  is 
biTugual). 

This  document  is  the  third  in  a  set  of  documents  that  provide  information  on  the  finite 
element  Model  for  the  Analysis  of  the  Vibrations  and  Acoustic  Radiation  of  Transducers 
(MAVART).  The  other  two  documents  in  the  set  are:  (1)  User’s  Manual  for  Program 
MAVART  and  (2)  Theoretical  Manual  for  Program  MAVART.  The  program  MAVART  has 
been  developed  for  the  Defence  Research  Establishment  Atlantic  (DREA)  under  research 
contracts  to  Canadian  industry  from  1976  to  the  present.  This  manual  provides  a  number  of 
sample  problems  to  aid  users  in  model  development  and  output  interpretation.  It 
demonstrates  virtually  aU  elements  and  features  of  MAVART.  It  also  provides  a  set  of 
problems  for  validating  modifications  and  new  installations. 


KtY  WUHUS,  utBUHII-'  I UHS  or  iDtN'l  ihlhHB  (teotinicaliy  meaninglul  terns  or  snort  phrases  that  aiaractenze  a  flocumara  and  could  be  helpful  in 
cataloguing  tha  document  They  should  be  selected  so  that  no  security  classification  is  required.  Identifiers,  such  as  equipment  modal  desfgiation,  trade  name,  military 
project  code  namejjeographic  locahon  may  also  be  included.  If  possible  keywords  should  be  selected  from  a  published  Ihesaunis.  e.g.  Thesaurus  of  Engineering  and 
Scientific  Terms  fItST)  and  that  thesaurus-identified.  If  it  not  possible  to  seled  indexing  terms  which  are  Unclassified,  the  classilicabon  of  each  shoulcfbe  indicated  as 
with  the  title). 

Underwater  Acoustics 

Transducer 

Design 

Analysis 

Finite  Element  Model 
Computer  Software 
Axisymmetric 
Piezoelectric 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 


INFORMATION  SCIENTIST'S  INITIALS 
INITIALES  DE  L'AGENT  D'lNFORMATION 
SCIENTIFIQUE  _ 


NO.  OF  COPIES 
NOMBRE  DE  COPIES 


COPY  NO. 
COPIE  N® 


AQUISITION  ROUTE 
FOURNI  PAR 


DATE 


DSIS  ACCESSION  NQ 
NUM^RO  DSIS 


1^1  National  Defense 
Defence  Rationale 

PLEASE  RETURN  THIS  DOCUMENT  PRI^RE  DE  RETOURNER  CE  DOCUMENT 
TO  THE  FOLLOWING  ADDRESS:  A  L'ADRESSE  SUIVANTE: 

DIRECTOR  DIRECTEUR 

SCIENTIFIC  INFORMATION  SERVICES  SERVICES  D'lNFORMATION  SCIENTIFIQUES 
NATIONAL  DEFENCE  QUARTIER  G^N^RAL 

HEADQUARTERS  DE  LA  DEFENSE  NATIONALE 

OTTAWA,  ONT.  -  CANADA  K1A  0K2  OTTAWA,  ONT.  -  CANADA  K1A  0K2 


